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Ylp-Chle  Lee*  and  Valter  Jaunzenla** 

Abstract 

The  objective  of  the  present  study  is  to  establish  thermodynamically 
validj  non-isothemal  stress  strain  relations  for  the  elastic -plastic 
rsuage,  and  to  obtain  solutions  of  typical  one -dimensional  problems  in¬ 
volving  unsteady  tanperatures . 

In  the  first  part  of  the  report  the  entropy  balance  equation  and 
the  expression  for  production  of  internal  entropy  are  considered.  The 
limitations  iii^)Osed  by  r^he  second  lav  of  thermodynamics  upon  plastic 
flow  rules  are  investigated,  £md  subsequently  a  set  of  non-lsothexmal 
plastic  stress  strain  relations  is  Introduced.  It  is  then  shewn  that 
these  relatlcxis  are  in  accord  vlth  theimodynamlc  irreversibility  when¬ 
ever  the  yield  criteria  of  Tresca  and  von  Mlses  are  used.  Moreover,  it 
is  foxmd  that  the  stress  strain  relations  of  von  Mises  are  a  special 
case  of  the  non-isotnezmal  stress  strain  relations  pre^osed  here.  An 
additional  problem  concerns  the  specification  of  elastic  unloading  from 
a  plastic  state,  and  in  the  presence  of  a  temperature  dependent  yield 
stress.  Precise  criteria  for  elastic  loading  and  repeated  plastic 
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flow  are  presented^  thus  coopletilng  the  characterisation  of  the  elastic - 
plastic  response. 

The  seccnd  part  of  the  report  is  conceded  with  an  appll^tlon  of 
the  general  theory  to  a  problem  for  the  infinite  half -space  constrained 
against  lateral  motion^  and  subjected  to  a  heat  p\tlse  unlfoxnly  applied 
over  Its  boundary.  In  the  analysis  of  this  problem  the  yield  criterion 
of  Tresca  is  tuaed;  the  treuislent  and  steady  stat6  solutions  for  both 
strain -hardening  and  perfectly  plastic  media,  having  either  a  constant 
or  a  ten^erature  dependent  yield  stress,  are  presented.  The  appearance 
of  elastic  and  plastic  regions  of  loading  and  unloading  Is  studied  In 
some  detail,  and  the  residual  stresses  and  deformations  are  correlated 
to  the  maximum  boundary  tenQ>erature . 

The  solution  of  the  half -space  problem  Is  then  extended  to  related 
problems  for  an  Infinite  plate  of  finite  thickness,  and  constrained  in 
the  lateral  direction.  It  is  shown  that  expressions  for  predicting  maxi¬ 
mum  residual  stresses  and  strains  can  be  obtained  directly  from  the  solu¬ 
tion  of  the  half -space  problem. 

Rumerlcsd  results  are  given  for  the  case  when  the  material  is  an 
aluminum  alloy.  In  addition,  the  maximum  boundary  tenperature  corres¬ 
ponding  to  the  rectirrence  of  plastic  flow  dxiring  a  second  heat  pulse 
Is  determined  Msisnlng  that  the  conditions  produced  by  the  first  pulse 
have  reached  a  steady  state. 

Althouc^  the  problem  selected  to  Illustrate  the  use  and  inplicaticns 
of  non-lsothermal  plastic  stress  strain  relaticnis  is  one-dimensional,  the 
method  of  approach  and  the  bMic  equations  may  be  used  in  the  analysis  of 
thermal  stress  problesis  in  two  or  three  dlmensicms. 
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PAKT  I.  TBESMOinnUNlCS  ABD  THBRfOFIASTICITY 


Oiapter  1.  Introduetloc  and  Review 
1-1  Introdnetlon 

Ibe  present  theory  of  plasticity  Is  applicable  only  to  problesu  in 
which  the  effects  of  tenqperature  on  the  plastic  stress  strain  relations 
aay  be  neglected.  In  particular^  the  yield  function  or  functions  which 
characterize  yielding  are  asstoned  to  be  independent  of  temperature. 

The  reason  for  the  neglect  of  non-isotheznal  probleuts  Is  the  relative 
coogtlexlty  of  the  plasticity  theory^  tdilch  would  be  further  increased 
by  the  inclusion  of  themal  effects. 

It  Is  well  known  that  the  yield  stress  of  most  engineering  aater- 
lals  decreases  with  Increase  In  tenperatxire  [1.1]*,  [1.2].  As  an  illus¬ 
tration,  we  note  that  for  Berrylliun  the  yield  stress  decreases  from  97 

3  5 

X  l(r  psl  at  room  teig>erattire  to  approximately  20  x  10  psi  at  1200  F, 

3 

whereas  for  Inconel  X  the  corresponding  decrease  is  from  112  x  10  psl 
to  99  X  10^  psl  [1.3].  If  the  Influence  of  temperature  <m  the  yield 
functl(»  and  on  the  stress  strain  relaticms  is  Included,  it  Is  possible 
that  the  stress  may  decrease  rather  than  increeuse  with  the  Increase  of 
total  strain.  Such  an  effect  would  have  no  counterpeurt  in  the  isother¬ 
mal  plasticity  theory. 

In  general,  plastic  flow  and  tenperature  change  occur  simultaneous¬ 
ly  and  influence  each  other.  The  search  for  the  possible  forms  of  non- 
isothermal  plastic  stress  strain  relations  Is  a  task  of  cmslderable 
theoreticeil  and  practical  value. 

*  luiibers  in  brackets  refer  to  the  Bibliography. 


1-2  Seope  of  tlia  Invgtlafctloii 


ms  report  is  esseatislly  divided  into  tve  perts,  Ifert  1  belag 
mlxUjr  eoDcemed  vlth  the  gmerel  theory  of  Bon>lsothezMl  plMtlelty. 
Qpeelfleally,  in  C!hapter  2  the  First  sad  Second  Levs  of  Ihemodpnenlesi 
together  vlth  the  concept  of  energy  conversionj  are  critically  discuss¬ 
ed  In  conjunction  vlth  plastic  flew.  She  restrlctloa  laposed  on  the 
non-lsotheznal  plastic  stress  strain  relations  hy  the  Second  Lav  is  in- 
vestigated>  and  the  conditions  of  loading  and  unloading  for  strain 
hardening  and  perfectly  plastic  naterlals  vlth  tenperature  dependent 
yield  stress  are  c«i8ldered  in  the  first  part  of  Chapter  lOn-lso- 
themal  plastic  stress  strain  relations  are  then  presented,  and  the  re¬ 
lations  associated  vlth  both  the  Nises  and  Tresca's  yield  criteria  are 
shown  to  satisfy  the  Second  lav  of  Thexmodynaaics.  Part  11  is  devoted 
to  the  application  of  the  non-lsothezaal  flow  rules  developed  in  Pasrt  1 
to  illustrative  problesu.  The  response  of  an  infinite  half  space,  con¬ 
strained  against  lateral  notion,  to  a  unlfoanly  applied  heat  pulse  over 
Its  boundary  Is  considered.  The  nedlun  of  the  half  space  is  assuned  to 
be  elastic,  linearly  stredn  hardening,  and  having  a  yield  stress  vhlch 
varies  llneenrly  vlth  teaperature.  Subsequently  the  cases  of  elastic, 
perfectly  plastic  medlim  with  a  constant  or  tes^erature  depaadrat  yield 
stress  are  discussed.  Flz»lly,  the  nethod  of  solution  is  applied  to  a 
plate  of  finite  thickness  in  order  to  obtain  an  estinate  on  the  dinen- 
sions  of  the  plastically  defomed  region,  since  in  the  half  space  prob- 
len  the  dlaenslon  of  length  rwsslns  arbitrary.  One  of  the  faces  of  the 
plate  is  assunsd  to  be  subjected  to  a  unlfom  heat  pulae,  vhersas  the 
other  face  is  either  insulated  or  naintained  at  aero  tanperature. 


.5 


luMrleal  eoagputatloiis  are  presentod  for  an  aluaima  alla9r»  aad  illat- 
trata  th«  tranda  la  tha  variation  of  tranalaat  and  raaldaal  atrataatf 
atralna  aad  dlsplaeauttitB.  Of  particular  Intarast  ara  the  dtMnslana 
of  plaatleally  daformed  regi«iB  as  datexnlnad  hy  the  aaplitadaa 

of  the  heat  pulse. 

1-5  Ravlair  of  Literature 

A  set  of  nan-lsothexnal  stress  strain  relations  for  plasticity  vas 
recently  proposed  by  Pnger  [1.^];  the  consistency  of  these  relations 
with  themodynsttlc  principles  has  not  been  investlflated,  hcwever  the 
vork  of  Prager  has  been  further  elaborated  by  Boley  [1.5]  and  Hat^dl 
[1.6]. 

Vakulenko  [!.?]>  [1.8]  and  Ziegler  [1.9]  have  both  attenpted  to 
foxBolate  the  non -Isothermal  plastic  stress  strain  relations  from  the 
point  of  vlev  of  Irreversible  thexmodynaalcs.  A  short  stpsary  of  the 
basic  concepts  of  this  subject  may  be  found  In  the  vrltings  of  Freuden- 
tfaal  [1.10],  Nac^idl  [1.11]  and  (Iregarlan  [1.12];  the  fundaaentads  of 
Irreversible  thezmodynsnics  axe  presented  in  the  vell-kncwn  texts  of 
Frlgoglne  [1.15]«  Callen  [1.14]  and  DeChroot  [I.I5].  Vakulenko  atnd 
Ziegler  start  vlth  the  assumption  that  the  genezaa.lzed  fozues,  e.g. 
stresses,  axe  Uneaurly  related  to  the  generalized  fluxes,  e.g.  plastic 
stradn  rates,  aklthou^  the  validity  of  this  assumption  Is  not  Imaediate- 
ly  obvious.  Both  authors  cocslder  the  plastic  power  as  the  only  qinan- 
tlty  associated  vlth  the  entropy  production  during  non-isothezmal  jtlas- 
tlc  deformation;  In  particular,  Ziegler  assumed  that  the  stress  aad  plas¬ 
tic  strain  rate  can  be  derived  from  tvo  potentiaLLs  the  sms  of  which  is 
equated  to  the  plastic  power.  For  a  strain  hardmilng  medium,  Vakulenko 


deriTed  a  flow  role  trm  a  plastic  potential  idiieh  was  foaad  to  l>e  a 
function  of  the  plastic  power,  the  free  energy  density  and  another 
stress  potential,  hut  laade  no  referwice  to  the  fom  which  the  flow  role 
would  asstone  during  unloading  and  neutral  change  of  state  ot  the  aediun, 
nor  was  there  any  mention  of  the  relation  between  a  yield  function  and 
the  flow  rule.  Although  hij^ily  interesting,  the  researches  of  Viakalen- 
ko  and  Ziegler  are  still  on  the  naujtoe  of  e^qploratoxy  attempts. 


J 


Chapter  2  thezmodyuaaics  and  Plastic  Oefoxnatlon 
2-1  Introdoetlon 

The  subject  of  thexnodjraaiiilcs  is  the  study  of  energy  and  entroipy 
dianges.  Sipeelflcally,  the  first  lav  of  thexnodynaalcs  Is  concexned 
vlth  energy  balance^  vhereas  the  second  lav  characterizes  the  Izrever- 
slblllty  of  defoxanitlons  through  the  Clausius -Oohen  lneq:aallty.  Any 
process,  reversible  or  Irreversible,  is  Invariably  governed  by  these 
tvo  lavs. 

Ve  shall  establish  In  vhat  foUovs  the  necessary  thexnodynaiaic 
concepts  and  develop  appropriate  foms  of  the  first  and  second  lavs  of 
thermodynamics.  For  convenience,  and  vlthout  loss  of  generality,  ve 
employ  rectangular  Cartesian  coordinates  throui^out. 

2-2  The  First  lav  of  Thermodynamics 

The  First  Lav  of  Thexmodynasilcs  expresses  the  condition  of  energy 
badanee  and,  specifically,  states  that  the  time -rate  of  change  of  the 
total  energy  in  a  body  is  equal  to  the  rate  of  stqpply  of  mechanical  and 
heat  energies.  We  denote  the  mechanical  power  of  external  forces  by  P 
and  ejqpress  the  flux  of  heat  into  the  body  D  as 

-/  “i  '*“  ■  -/  ’i.i " 

B  D 

vhere  q  Is  the  heat  flux  vector,  n  Is  the  outer  noxml  of  the 
boundary  B  of  the  region  D,  and  do,  dr  are  surface  and  voluae 
elements  respectively. 

The  First  lanr  may  then  be  stated  as  follows. 


dT  ■  K 


$ 


-/  *1,1 


where  K  Is  the  material  time-rate  of  the  klaetic  energy  K,  and 
Is  the  material  time-rate  of  the  total  internal  energy. 

Since 


K .  i  y'  p^i  ,j  dt  , 


where  f>  is  the  mass  density^  and  the  particle  velocity,  it 
follows  that 


K 


/  “i  ’i 


dT 


(b) 


Here  a^  denotes  the  particle  acceleration  defined  by 


The  mechanical  power  F  is  given  by 


ijTjdd  *  y'  pfjvj 


d  T 


(c) 


where  and  f^  are  the  coaoponents  of  the  STirface  tractions  T 

and  body  force  f,  respectively.  The  stress  tensor,  and  the 

stress  vector  T^  are  related  by  the  equation 


•^i  ■  ’ij  “j  ' 


^diich,  when  substituted  into  (c)  yields 

^  ■/  “  ’ji,i  *  “i>  ’'i  *  ^1  'i.j'  * 


(4) 


•7 


after  ualag  the  dlTergenee  theoren  for  the  lerface  Integral.  Ylrtue 
of  the  eq^tlon  of  aatlan> 

O'!-  I*!  <•> 

(d)  reduces  to 

i’  J  ’ij <') 

D 

fOnaaUjr,  ve  deeonpose  the  velocity  gradients  v.  .  into  syaaetric  and 
amtisyieaetrie  parts, 


+  V 


ij 


vhere, 


iJ 


'ij 


I  ^^ij  ■ 


(g) 


are  referred  to  as  defcmation  rate  and  vorticity,  respectively. 
From  (h)  and  (f)  now  follows 

D 


Since  the  stress  tensor  is  synnetric,  ve  have  that 


T  V  s  0 

iJ  iJ 


Therefore,  letting 


A  ■  J'  p  X  if 

D 


(1) 
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vhere  X  1b  tbe  rate  of  Internal  energy  density,  Ye  arrive  at  the  foxn 

[2.11,  t2.^,  [2.51,  [a.*] 

PI  *  %,1  -  ■>„  (2'2-2) 

Of  the  lav  of  conservation  of  energy.  It  Is  pladslble  to  write 
X=U+j^  (J) 

« 

idiere  U  Is  the  time  rate  of  a  function  U  of  elastic  strain 
and  the  entropy  density  s,  whereeis  ^  represents  the  time  rate  of 
internal  energy  due  to  plastic  deformation,  and  Is  not  the  time  rate  of 
a  state  function  since  pleistic  defonnatlcm  is  a  path  dependent  process. 
Then 


»  ”  *  O  I' ‘  S.,!  ■ 

Eq,uatlon  (2.2.3)  Is  a  representation  of  the  principle  of  bedance  of 
mechanical  and  thermal  energies  in  a  fom  that  is  most  c(xivenlent  for 
subsequent  investigation. 


2-3  Heat  Conduction  Equation  and  the  Entropy  Balance  Bamtlon 

Ve  assume  that  the  energy  U  in  equation  (j)  of  the  preceding  sec¬ 
tion  is  a  function  of  the  elastic  strain  and  the  entropy  s,  thus 


U  -  U  ,  s) 


(2.3.1) 


Ve  also  adopt  the  customary  form  of  the  equations  of  stats  [2.3], 

[2.6], 

(2.3.2) 


ij 


OS 


In  the  ease  of  wall  defoiaatlons  and  taqpwratore  ehaaoesi  the 
funetloa  U  any  be  approximated  by  the  first  teias  of  a  Saylor's  ex- 
pansloi  about  some  reference  state.  Taking  the  values  of  U, 

T  and  8  at  the  reference  state  to  be  zero,  ve  obtadu 


U  •  (As  +  2b  sa^  +  j-J"  +  2|ji  ^ji^ 


(2.3.3) 


Sy  (2.3.2)  and  (2.3.3),  it  foUcws  that 


T  *  (a  8  +  b/^  )/p 


’ij  •  p  “  'll  *  ^ 


equivalently. 


8  =  (Tp  -  bA^)/A 


^ij  ®ij  "  Ap^  ^  A  ®ij 


(2.3.4) 


where  A,  b,  7  in  (2.3.3)  and  (2.3.4)  are  known  material  properties 


and  is  equal  to  We  write 

„  -g,. . 


(a) 


-  J  -  (31  +  2p)  a 


(!») 


where  X  and  p  are  Isuse  constants  and  oc  is  the  coefficient  of 
thexnal  expansion. 


(2.3.5) 


Be^Ulngttat  djj  1.  th.  totol  d.(«-tloa  mt.,  «  !>.« 

*1J  ■  ^IJ  *  ‘ij 

vhere  6. ,  and  8  .  are  the  elastic  axid  plastic  defoznation  rates, 

XJ  IJ 

respectively.  Substitution  of  (2.3.2)  and  (2.3.5)  in  (2.2.3)  th«i 
yields: 

ELimioating  a  in  (2.3.6)  by  the  first  equation  of  (2.3.4)  and  replac- 
ing  ^  1  ^  ii  the  use  of  the  Fourier's  Lenr  of  Heat  Con¬ 

duction,  ve  arrive  at 

e 

p  G(T)  T  -  k  =  3^  +  (  s^j  -  p  ft)  (2.3.7) 


vhere  G(t}  »  T  p/A.  In  particular  for  small  temperature  changes  and 
defomation  ve  have  that 


2  m  ^ 

Po  ^o^ 


/v>  -  k  T 
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Ve  take 


V 


(c) 


vhere  p^  arre  the  Initial  values  of  the  mass  density  and  teBqpezmtore 
respectively,  c  is  the  specific  heat  per  unit  volume  and  k  is  the 
thermal  diffusivity.  Ve  new  havetUswe  conditions,  (a),  (b)  and  (e)  for 


u 


explicitly  det<pnlnlng  k,  t  end  j. 

The  ejqpressloi  given  by  eq^tlon  (2.3.7)  !■  knovn  ee  the  heat  coa- 
duetioa  equation.  The  first  tezn  on  the  rij^t  head  side  of  (2.3«7)  vsy 
he  Intexpreted  as  a  heat  source  due  to  the  rate  of  elastic  volune  defor 
nation,  vhereas  the  term  in  the  parentheses  represents  another  heat 
source  accounting  for  the  dissipation  of  the  inelastic  mechanical  ener¬ 
gy  into  heat  during  plastic  defoznatlon.  Since  plastic  volume  change 
is  usually  aussumed  to  he  zero,  the  tvo  heat  sources  represent  two  dis¬ 
tinctly  different  phenomenons.  The  tern  Tb^z/A  represents  a  second- 
order  effect  which  becomes  noticeable  only  when  the  duration  of  motion 
is  very  long.  Neglecting  this  term,  we  write 
« 

p  0(T)  T  -  k  s^j  -  p  JJT)  (2.3.8) 

It  should  he  noted  that  in  the  dissipation  of  inelastic  mechsmlcal  en¬ 
ergy,  not  all  of  the  plastic  power  s^^j  Is  converted  into  heat; 
only  the  part  in  excess  of  the  work  fi/jl  required  to  increase  the  inter 
nsl  energy  in  conjunction  with  plastic  deformation  is  dissipated.  The 
function  p)/  may  he  considered  as  the  rate  at  which  energy  is  absorbed 
by  the  medium  to  change  its  internal  structure  during  plastic  deforma¬ 
tion. 

Dividing  (2.3.6)  by  T,  and  making  use  of  identity 

we  obtain  the  equivalent  form: 

P  i  +  ^  ^  (’U  '14  - 


P^)  (2.3.10) 
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known  as  the  entropy  balance  equation. 

Let  us  new  Introduce  an  extezml  entropy  flux  S  by 

K  •/  -  ■/  h 

D  B 


D 


^ere  is  the  corresponding  density.  Writing 


(2.3.11) 


s  -  s  +  s. 
e  1 


(2.3.12) 


where  s^  is  defined  as  the  Internal  entropy  production.  We  note  that 

by  (2.3.10)  and  (2.3.11), 


P  s 


1 


P  8  +  ( 


(2.3.13) 


The  Second  Law  of  Thermodynsmlcs  In  the  fonn  of  the  Claxislus-Duhem  in¬ 
equality,  requires  that  for  any  process,  the  Internal  entropy  produc¬ 
tion  must  be  non-negative  (2.7],  [2.8],  [2.9]. 

P  8  +  (■-^),i  ^  0  (2.3.14) 


This  implies  that  the  rle^t  hand  side  of  (2.3.10)  must  be  greater 
than  zero.  Since  each  term  on  the  rl^t  hand  side  of  (2.3.10)  repre¬ 
sents  en  independent  phenomenon,  we  require  that  these  terms  be  scpsur- 
ately  non-negative,  thus 


-  % 


>  0 


(2.3.15) 


1? 


and 

By  Fourier  law  of  heat  conduction 

^  ^  ^,11 

Equation  (2.3>ll)  therefOK  requires  that 

k>  0 

The  Inequality  represented  by  {2.^.l6)  must  be  satisfied  In  any 
plastic  process.  It  may  be  further  explored  by  examining  the  heat  con¬ 
duction  equation^  (2. 3 >8)  for  two  extreme  cases**  If  all  the  plastic 
pover  Is  absorbed  by  the  material,  there  Is  no  dissipation  of  mechani¬ 
cal  energy  Into  heat.  This  vould  Isply  that  the  heat  source  in  (2.3*6) 
slnqply  vanishes,  hence 

‘ij  (2-5-17) 

On  the  other  hsuid.  If  no  energy  Is  being  absorbed  by  the  material  dur¬ 
ing  plastic  defoimatlon,  all  the  plastic  porer  Is  converted  Into  heat, 
hence 


p  -  0  (2*3*18) 

In  reality  these  tiro  extreme  cases  are  not  likely  to  obtain  [2*10]* 
It  has  been  demonstrated  by  Hort  [2.11],  [2.12],  Sato  [2.13],  Taylor, 
Farren  and  ()ainney  [2*14],  [2.13],  [2.l6]  using  calorimetric  measure¬ 
ments  of  tension,  torsion  and  ccsqiresslon  specimens  that  not  all  the 


plaatlc  pcwer  is  converted  Into  heat.  For  exsople,  Tagrlor  ai^  Farr«i 
[2.17]  found  In  their  e^eriaenta  that  for  steel  of  the  plastic 

vork  was  dissipated  in  the  form  of  heatj  hence  only  vas  absorbed 

by  the  material.  The  percentage  of  the  plastic  vork  converted  into 
heat  hM  been  found  to  range  from  9O.5  to  92  for  aluminum  and  from  92 
to  93  for  single  aluminum  crystal. 

We  may  postulate  that 

^  I*  ■  *.  ’ij  'ij 

idiere  is  a  property  of  the  material  which  can  be  experimentally 
determined  for  a  particuleur  material  euid  assumes  values  between  zero 
and  cne,  l.e.> 

1>A^>0  (2.3.20) 

The  cases  of  A  being  equal  to  zero  and  one  vinild  correspond  to 
e 

the  two  extreme  cases  represented  by  (2.3.17)  and  (2.3.I8).  Combining 
the  relations  (2.3.I6)  and  (2.3.19),  we  obtain 

I  'ij’lj 

Since  T  is  always  greater  than  zero,  it  follows  from  (2.3.20), 
that  the  requlrenent  (2.3.21)  is  eqpiivalent  to 

hi  •ij>  ° 

We  have  therefore  shown  that,  on  the  basis  of  the  assusqption 
(2.3.19),  any  flow  rule  which  satisfies  relatlcm  (2.3.22),  will  also 
satisfy  the  relation  (2.3.16)  obtsdned  from  the  Second  law  of  thermo- 
dynaalcs. 


Chapter  3  HoQ-Ieotheznal  Plastic  Stress  Strain  Belations 

3-1  Yield  Function  and  the  Criterion  of  Iioad1ng«  Ohloading  and  Meetral 

Change  of  State 

Itae  internal  forces  of  a  continuous  atedixm  are  described  by  a 
stress  tensor  Representing  the  state  of  stress  by  a  point  in 

six -dimensional  space  of  a  rnediun  is  custonmurily  said  to  be 

elastic  if  the  state  point  lies  vithin  a  certain  e«ivex  dcnain  contain¬ 
ing  the  origin  of  the  stress  space.  This  domain  is  referred  to  as  the 
elastic  domain^  and  ve  say  that  yielding  occurs  aus  soon  as  the  state 
point  reaches  the  bovindary  of  the  elastic  domain  characterized  by 


r  (*ij,  T,  t)  =  0 


(5.1.1) 


Here  the  yield  function  F  is  assumed  to  be  a  function  of  stress  x 
temperature  T,  emd  of  a  strain  hardening  parameter  (  defined  by 


ij 


} 


t 

*  ■/  <5.1.2) 

o 

One  of  the  basic  assusiptions  of  pleustlclty  theory  is  that  plastic 
flow  is  Independent  of  hydrostatic  pressure.  Therefore,  ve  shall  re¬ 
place  in  (3*l<l)  the  stress  x. .  by  the  devlatorlc  stress  v . .  defln- 

XJ  IJ 

ed  by 


'^IJ  “  "^ij  "  3  ®ij 
The  equation  (3>l.l)  new  becosies 


(3.1.3) 


T,  5  )  -  0 


(3.1.k) 


11$ 


and  va  lat 


F  1)  -  f  (v^j)  -  H  (T,  t) 


(3.1.5) 


vhere  }  Is  a  function  of  deviatorlc  stress  alone  and  H  Is  a 

IJ 

function  of  teoqperature  T  and  the  strain  hardening  paraawter  (. 
Since  the  aedlrDS  is  sussumed  to  he  hosiogeneous  and  Isotroiplc,  the  func- 
tl(«i  f  de9;>ends  on  throu^  the  Invariants 

^2  ^"^ip  “  I  "^IJ  ^IJ  “  I  ”^2  ^3  ^ 


(3.1.6) 


^3  ^^IJ^  "  5  ^ij  ^ki  “  5  ^  ^ 


Both  functions  f  and  H  are  further  assumed  to  he  continuous  and 

differentiahle.  Ve  shall  adopt  as  a  basic  postxilate  the  condition  that 

H  Increases  monotanlcally  with  respect  to  the  strain  hardening  para¬ 
s'’ 

meter  l.e. 

>  0  (3.1.7) 

Most  engineering  materials  exhibit  the  characteristic  that  the 
yield  stress  decreases  with  increasing  tenperature  [3.2],  [3.3].  Ve 
ohseirve  from  (3.1.^)  and  (3.1.5)  that  during  plastic  flow  the  boundary 
of  the  elastic  domain  expsmds  as  a  result  of  strain  hardening;  and  that 
this  boundary  may  sJso  ccmtract  or  espand  depending  on  the  rise  or  fall 
of  the  temperature.  If  the  Influence  of  the  temperatu3re  on  the  yield 
stress  donlnates,  it  is  probable  that  with  increasing  temperature  the 


*A  similar  postxilate  is  used  by  Hill  [3.1]  for  the  ease  when  the  strain 
hsxdenlng  parameter  is  equal  to  the  total  plastic  work. 
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elMtlc  dcnaln  nay  still  be  shrljaklng,  even  If  there  Is  strain  baxdcn- 
Ing.  Thus  the  possibility  of  a  decreatse  In  stress  aecoqpanled  by  an 
Increase  In  strain  and  teo^erature  In  the  nedlun  cannot  be  ruled  out. 

In  Isotheznial  plasticity,  a  medium  that  exhibits  this  phenoaienon  of  de¬ 
crease  In  stress  vlth  Increase  In  strain  Is  considered  sis  unstable 
since  It  behaves  like  am  energy  source.  Furthermore,  In  Isothermeil 
plasticity  the  decrease  In  stress  automatically  Inqplles  tmloadlng, 
which  Is  not  always  the  case  In  non-lsothexnHil  plasticity.  It  Is, 
therefore,  necessary  to  examine  the  Influence  of  temperature  on  the 
yield  functlcm  with  care,  amd  determine  the  conditions  under  which  the 
medium  undergoes  loading,  xmloeidlng  or  neutral  change  of  state. 

As  Is  customary  In  plasticity  theory,  we  set  the  yield  function 
equeil  to  zero  during  plaustlc  deformation.  If  the  stress  point  moves 
from  one  plastic  state  to  smother  In  stress  spsme,  the  plastic  strain 
Increases,  and  we  have  the  conditions 


F.O  ,  F.O  ,  0  ,  S>0 


The  vsmlshlng  of  F  implies 


JL-  V  ^  T  +  ^  »  n 
^  Vij  +  ^  T  +  S  =  0 


(5.1.8) 


By  (5.1.5)f  noting  that 


dF  Bo 

ST  -  -  5? 


(5.1.9) 


l8 


The  eqiuAtlOQ  (3.1.8)  nov  l)ecomes 


ar 


ij 


.  ar 


dH 


(3.1.10) 


Sines  (  is  a  positive  definite  quantity,  as  Is  eleeur  from  the  deflnl> 
tlon  (3.1.2),  and  recall Ing  the  hypothesis  (3.1.7)«  It  foUoirs  that 


dH 

¥ 


i  >0 


(3.1.11) 


The  criterion  of  loading  is  thus  es^ressed  by 


and 


p.o  ,  ;  =  o  ,  i>o 

^  *  S  (5-1.12) 


During  unloading,  the  plastic  strain  rates  become  zero  and  the  stress 
point  moves  from  the  yield  sxirfacc  bach  into  the  elastic  dosaiin.  The 
criterion  of  unloading  is  ejcpressed  by 


r<o  or  7.0  ,  7<0  ,  .JJ-O  ,  1-0 


and 


^  T  <  0 

5^  11  s 


(3.1.13) 


The  neutral  change  of  state  Is  then  defined  by 


and 


P«0,  F-0,  6"® 


ar  *  ^  fip  * 

5^  ^ij  ^  ^ 


T  -  0 


(3.1.14) 
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For  B  perfectly  plastic  materlad,  the  yield  function  F  is  a 
function  of  stress  v. .  and  teitQ>erature  T  alone,  i.e., 


F(v^j  ,  T)  -  f(v^j)  -  H(f) 


(3.1.15) 


the  criterion  of  loading,  unloading  and  neutrsLL  change  of  state  can  he 
deduced  from  (3.1.15),  (3.1.10),  (3.1.12),  (3.1.13)  and  (3.1.1^)  hy 
setting  equal  to  zero.  The  results  are  as  follows: 

Loading 


and 


F=0,  F  =  0  ,  8^j$0  ,  6>0 


dF 


(3.1.16) 


5——  V .  .  +  ^  T  =  0 

5;^  ij  s 


Utaloadlng 


F«=0  or  F  =  0  ,  F<0  ,  s^j  =  0  ,  I  -  0 


and 


^  i  T<o 

^  ±j 


(3.1.17) 


Reutral  change  of  state 

F=0  ,  F=0  ,  s 


ij 


S  -  0 


dF  *  dF  *  « 
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(3.1.18) 


3-2  Development  of  Koa-IsothentiaJ.  Plastic  Stress  Strain  Relmtiop* 

In  this  section  both  the  strain  haLrdenlng  and  perfectly  plastic 


media  vlth  tan^erature  d^endent  yield  stress  are  considered.  Since  the 
elastic  response  Is  characterized  by  Hooke's  lav  vlth  an  added  tempera¬ 
ture  term,  attention  la  directed  mainly  to  the  non -Isothermal  plastic 
stress  strain  relations.  It  Is  assimed  that  an  Infinitesimal  change  of 
the  plastic  strain  r. .  Is  caused  by  Infinitesimal  change  In  the  devla- 
torlc  stress  t.  .  and  the  tesqperature  T.  Furthermore,  the  relations 
betveen  the  rates  of  plastic  strain,  stress  and  temperature  must  be  homo¬ 
geneous  of  order  one  so  that  any  results  obtained  vould  be  Indepedent  of 
the  time  scale  en^loyed.  We  vrlte 


®1J  “  \jki  "^kl  ®1J 


(5.2.1) 


vhere  emd  depend  on  T  and  r^j.  Incorporating 

the  conditions  (3.1.12),  (3.1.13),  (3.1.1^)  for  loading,  unloading  and 
neutral  change  of  state  Into  (3.2.1),  the  coefficients  emd 

teike  on  the  form 


A  =  C  ^ 
*ljkl  ^IJ 


B 


r  ^ 
ij  “  ‘^ij  Sr 


(3.2.2) 


and  the  flov  lav  (3.2.1)  becomes 


‘ij  ■  '  0  (5-  y  ■  ^  - 


If  F  ^  0  ,  F  <  0 
0 


Vy.,  +  ^  T)  If  F  -  0  ,  F 


Cxa/33 


In  t1«v  of  (3«l<I>2)j  tho  Quantity  in  tlw  paronthoiit  in  (3*8<9)  la  al- 
naya  poaltlva  daflnlta  dturlng  plaatie  flov.  fha  plaatle  atraaa  atxmln 
ralation  anat  further  satisfy  tvo  nseossazy  eondltlona  aaaaly, 

tha  eondltlon  of  plastic  ineosipresslhillty, 

8  -  0 
11 

and  ths  oondltloa  that  the  direction  of  principal  axes  of  plastic  strain 
Increaent  coincides  vlth  the  direction  of  the  principal  stress  axies« 
since  the  nedlum  is  asstimed  to  be  isotropic.  These  tvo  conditions  are 
satisfied  by  choosing 


vhere  the  function  g  my  be  considered  as  a  plastic  potential,  and  de« 
pends  on  the  invariants  and  Sy  vhereas  D  is  a  function  of 
T  and  |.  The  most  convenient  choice  of  the  function  g  is  to  Identi¬ 
fy  it  vlth  the  function  f  introduced  in  (3.1.5)*  Rspiaclng  g  by  f 
In  (3. 2. If)  ve  obtain 


C.,  -  i  ^  ^  )  (5-2-5) 


'IJ  D 


IJ 


'2  r'lj 


'3  ''•ij 


Fron  the  definition  (3*1*6}  of  and  it  Is  found  that 


^2  ^ 


ij  W'lj 


9 


(3.2.6) 
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From  (J.2.6)  and  (3*2. ?)  ve  readily  deduce  that  ■  0,  and  thua  the 
condition  of  plaetic  ineonpressihility  is  net.  In  terns  of  the  principal 
rates  of  stress  and  strain,  the  relation  (3*2.3)  nay  be  written  as 


*i  “  D  ^  ^  \  ^  M  ° 

where  s^,  and  are  principal  plastic  strain  rates,  principal 
total  and  devlatorlc  stresses,  respectively.  The  cosnonent  df/dr^ 
represents  the  projection  of  the  gradient  of  the  function  f  in  the 
direction  of  the  prlnciped  stress 

It  remains  to  determine  the  function  D.  Recalling  that  during 

»  • 

plastic  flair  both  F  and  F  are  equal  to  zero  and  noting  that  (  is 

1/2 

equal  to  (s^j  s^^)  ,  it  is  found  after  substituting  (3.2. 3)>  (3*2.5) 

into  (3*1.10)  that 


We  recall  that  wm  postulated  to  be  positive  definite  during 

plastic  flew,  and  observe  that 

1/2 


is  positive  definite  if  the  positive  root  of  the  quantity  in  the  parM- 
thesis  is  taken.  Therefore  the  quantity  D  is  positive  definite. 

Coabinins  (3*2.3),  (3*2.5)  and  (3*2.7);  the  non-isothexnsil  plastic 


flow  law  new  becoows 


*±i 


f or  F  ^  0,  F  4  0 


dH  /  ^ ^ _ N 


w 


kl 


(3.2.8) 


The  foxnulA  (3.2.8)  Is  a  general  flov  rule  for  strain  hardening  media 
vlth  tengperature  dependent  yield  stress,  and  can  be  used  to  determine 
the  plastic  strain  rates  once  the  function  F  Is  knovn. 

The  tvo  most  commonly  used  yield  criteria  are  due  to  Von  Mlses  and 
Tresca  respectively.  The  yield  function  defined  by  the  Von  Mlses  yield 
criterion  is  given  by  [5.5]>  [3.6] 


F  XT  f  -  H 

=  Jg  -  (T,5)  (3.2.9) 

where 

f  =  Jg  ,  H  =  (T,|) 


and  K  Is  the  yield  stress  In  slaiple  shear.  The  yield  function  repre¬ 
sented  by  (3.2.9)  Is  said  to  be  regular  because  a  unique  normal  can  be 
defined  at  every  point  on  the  yield  surfhce.  If  (3.2.9)  Is  substituted 
Into  (3.2.8),  the  flow  rule  associated  with  the  Von  Mlses  yield  criter¬ 
ion  then  tadces  on  the  form 


®ij 


dH 


(v^ 
mn  mn 


for  F  <  0,  F  <■  0 

^ - 375“  (  ^  ^  for  F  -  0,  F  -  0 


) 


(3.2.10) 


o 


She  flow  rule  for  a  perfectly  plauitlc  medium  vlth  teagperature  dependent 
yield  strees  can  be  established  by  considering  aigain  the  condition  of 
the  coincidence  of  the  direction  of  principed.  axes  of  plastic  strain 
rates  vlth  the  principal  stress  euces  during  plutlc  flow,  and  the  con¬ 
dition  of  plastic  InccBpresslblllty.  These  conditicms  suggest  the  flow 
rule 

0  for  F  -  0  ,  F  ^0 

Sf  ,  (5-a-ii) 

for  F  ■  0  ,  F  ■  0 

vhere  X  is  a  function  of  plastic  strain,  stress  and  tenqperature  and 
their  rates.  It  is  not  possible  to  detexmine  \  from  the  yield  fono- 
tlon  Itself  because  the  latter  Is  Independent  of  plastic  strain.  If  ve 
consider  the  yield  criterion  of  Von  Mlses, 


F  =  f  -  H 


-  K^(T) 


(3.2.12) 


then  X  may  be  obtained  by  multiplying  (3.2.11)  by  and  sunning. 

In  view  of  the  definition  (3.1.6)  of  Jg,  ve  find  that 


X  = 


(3.2.13) 


3“3  Hon-Isothemal  Plastic  Stress  Strain  Relaticms  for  Sing»T*^'»*  Yield 
Functions 


The  yield  function  defined  by  the  Von  Nlses  yield  criterion  is  a 


iS 


xtgultr  yield  foaetioa  representing  a  yield  eurfaes  without  eomers.  Ve 

My  also  oharaoterlse  yielding  hy  a  set  of  h  yield  foaetlons 

C^lJ'  where  7  ■  l^&>...a.  these  yield  fanotlOns  define  a 

yield  surface  in  the  stress  space  which  My  not  hawe  a  unlq;as  at 

every  point,  and  are  referred  to  as  singular  yield  fonctloas.  Ve  adqpt 

the  convention  that  plastic  flow  begins  whenever  any  (me  of  the  fune> 

(7) 

tlons  T  beoones  zero,  and  that  the  aedlua  is  elastic  if  all  fane- 
(7) 

ti(ms  F  are  less  than  zero.  The  plastic  strain  rate  is  then  taken 
to  be  the  sun  of  the  constituents  [3.71,  [3.8],  [3.9],  tr>*iw 


•ij 


i  s 

A  « 


(3.3.1) 


For  a  strain  hardening  meditsn  with  a  tenperature  dependent  yield 
stress,  the  plastic  strain  rates  are  defined  by 


s 

IJ 

-  0 , 

iW  - 

0  ,  F^’'^ 

• 

8,(7) 

• 

m  ^  0 

IJ 

1  ^  U 

s 

■ij 

t=~o 

^iJ  " 

^(7) 

e 

T>  0 

(3.3.2) 


;(r) 


where  end  are  the  values  of  s^.  and  (  oorrespooding 

(y) 

to  the  function  F'  '«  The  most  consaonly  used  singular  yield  functions 
due  to  Tfesca,  and  are  defined  by  [3.10],  [3.11] 
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y(l)  ^  ^(1)  ^ 

,(2)  ,  ,(2)  _ 

p(5)  „  f(3)  _  h(T,|) 

vhere 

,(i)  .  1-^1  •  ’2i  ,  Pi  - 

^22 

,(2)  .  P2  -  V  .  1^2-V 

2  2 

p)  .  P;  -  M  .  P3  - 

^  2  2 

amd  are  the  prlncipcLL  mlues  of  the  devlatoric  and  total 

stress  tensors.  The  function  H  is  identified  with  the  yield  stress; 
for  a  perfectly  plastic  medium  H  is  a  function  of  tenperatura  T  aLLcne. 
The  plastic  strain  rates  to  (5.3,1)  (3.5.3)  |:^1»lpa 

ed  from  (3.2.8)  by  replacing  the  functions  f,  P  and  H  by  the  corres¬ 
ponding  functions  f^^^,  P^^^  and  H  associated  with  the  singular 

yield  surface. 

Ve  note  that,  for  Tresca's  yield  criterion  not  more  than  two  of  the 
three  singular  yield  functions  may  vanish  [3.13]  ,  the  third  one  would 
be  equal  to  -E.  In  this  case  the  stress  point  is  at  the  comer  of  the 
yield  surface  defined  by  (3*5.5)  in  the  space  of  principal  stresses.  Ve 
renarh  that  Tresca's  yield  criterion  is  frequently  used  when  the  direc¬ 
tion  and  sense  of  the  priucipcd.  stresses  are  known  from  cooslderatlon  of 
syaawtry. 

Prom  (5.2.11)  follows  that,  for  a  perfTOtly  plastic  aeditsi,  the 


(3.3.5) 


(3.3.i^) 


flow  rule  ueoeiated  with  the  Tresoa's  yield  eriterion  le  ej^reaMd  Iqr 


(3.3.5) 


where 


-  0  if  0  or  -  0  ,  P^’'^< 

ap<y>  • 


if  P^’'^ 


ar  •  3f  :  . 

s;^  \  ^  ^  ®  ° 


(3.3.6) 


She  quantities  X  are  functions  of  stress,  teBg>erattDre,  strain  and 
their  rates.  !Ehese  coefficients  may  he  determined  in  a  manner  similar 
to  that  used  in  the  last  section,  althouf^  for  particular  problens  other 
methods  may  be  more  suitable  (cf.  chapter  ^). 


3-4  The 


c  Investiflatlon  of  the  Hcsr  Rule 


The  conditions  iiqposed  on  fleer  rules  by  the  principle  of  thermo* 
dynamics  were  derived  in  the  preceding  chapter.  Specifically,  ve  found 
that  nan<iegatlve  entropy  production  requires  that 


^IJ  ®ij^  °  (2.3.22)B 

be  valid.  Ve  Intend  to  sheer  in  this  section  that  the, general  non-lso- 
themal  fleer  rule  given  by  (3.2.8)  satisfies  (2.3.22),  and,  thersfoire. 


.89 


the  second  lav  of  thexnodynanles,  idxen  used  in  conjunction  vlth  either 
the  Von  Mlses  or  Tresca's  yield  criterion. 

Ve  recall  the  flew  rule  (j.2.10)  associated  vlth  yield  criterion  of 
Von  Mlses,  and,  multiplying  It  by  and  simnalng,  obtain  the  result 


''iJ  ®1J 


^IJ  ^IJ 


dE 


^5.4.1) 


In  order  to  verify  (2.3.22),  we  first  recall 


(3.1.12)R 


and  that  dH/d(  was  also  assumed  to  be  positive  (cf.  relation  (3>1«7)’)* 

1/2 

Then,  since  ^mn  ^  clearly  positive,  the  only  quantity  that 

needs  to  be  examined  is  t  .  .  v.  . . 


Writing 


T 


ij 


IJ 


+ 


1 

5 


(5.1.5)R 


ve  find  t^t 


T.  .  V.  »  V.  ,  V,  .  +  -T  B, .  V.  .  T. . 

iJ  iJ  iJ  iJ  5  iJ  iJ  kk 


’'ij  ’ijs  ° 


lAere  ve  used  the  fact  that,  be  definition 


'11 


(3.4.2) 


0 


•aj 


i.e.>  the  first  InT&riant  of  the  devlatorie  stress  tensor  -vanishes  ld«i- 
tloally.  fhe  result  (3.4.2)  thus  cGsqpletes  the  proof  of  (2.3.22)  tor 
the  ease  -when  the  yield  criterion  of  Von  Hises  is  adopted. 

In  order  to  verify  (2.3.22)  for  -the  flow  rule  (3*2.8)  used  in  con¬ 
junction  with  the  Tresca's  yield  criterion,  ve  consider  a  stress  point 
at  a  comer  of  the  yield  surface.  Qie  case  when  only  one  singuJar 
function  vanishes  is  then  included  eis  a  special  ease.  Let  us  assvie 
that  in  the  yield  comer 


V  >  V 

1  2 


(3*4.3) 


and  that  the  singular  yield  functions  are  now  r^resented  by 


<1)  .  -H.o 

<2)  .  Xp.  .H.o 
(3)  .  hpA  .H.o 


(5.*.  4) 


In  addition,  ve  let 

p(l)  =  0  ,  .  0  ,  <  0 


(3*4*5) 


Noting  that 


■^1  "  “^2  "  “^1  ■  "^2  '  "^l  "  "^3  “  ‘'l  ■  *^3 


the  assumption  (3*4*3)  iimplies 


^1  -  *^2^  0  '  ■'l  -  0 


(3*4*6) 
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Ihe  first  tvo  eqiostloiis  of  esp  be  vrlttea  ms 


,(i)  .  .H 


.(a)  ''i  •  %  ^  ■  ’5  . 

f'  '  ■  . A  •  H 


Since  the  function  H  is  the  sane  for  all  three  singular  yield 
functi(Xi8,  this  ingplies  that 


Vg-Vj 


^2  "  “^5 


(5A.8) 


where  f®  and  P®  are  the  values  of  the  functions  f^^\  f^^^  and 

P^^^,  P^^^  at  the  yield  corner*  respectively. 


Multiplying  (5.3*l)  by  and  noting  that 


^ij  “ij  “  ^  “k 


ve  obtain 


^k«k 


“k 

7-1 


(7) 


(3.J^.9) 


Since*  by  (3.2.8)* 


(7) 


ae  / 


T)  (3.4.10) 
'”n 
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the  evibetltutlon  of  (3*4. 10)  Into  (3>^>9)  vltb  (3«^*3}«  (3«^>7)  *uad 
(3<h.8)  yield* 


Vl  ^  V2  V3 


'^1  <ir* 


— TT^  ^ 

iE  ,it‘  it‘  ” 


2 

ap® 


)] 


T)  (3.4.11) 


a  a 


) 


The  flov  rule  associated  with  the  Tresca's  yield  criterion  is  com> 
patlhle  vlth  the  seccsid  lav  of  thermodynamics  if  the  ri£^t  hand  side  of 
(3.4.11)  can  be  shown  to  be  greater  than  zero.  The  quantities  (dai^0> 
(bf^/br^  bf^/br^)  ^  and  (dP^/av^^  +  aP®/aP  T)  sure  adi  positive 

definite  during  plastic  flow.  The  only  quantity  that  needs  to  be  exam¬ 
ined  is  the  sum  of  the  terms  Inside  the  square  bracket,  fiy  (3*4.7)  <uad 
(3.4.8),  we  obtain 

df^^^  Sf^^^  1  df^^^  1  df^^^  1 

^  “  5?;^*  2  '  ^ - 2  »  - 2 


(3.4.12) 


Hence,  in  view  of  (3.4.6),  ve  have  that 


sa 


[  ( 


g -  +  S -  )  +  ^2  (  - 

fjT  dT  2  «. 


^rl  ('i  -  'a)  ♦  r  (’1  -  'slj  ^ 


Thus  It  has  been  shovn  that  the  flew  rule  used  on  conjunction  with 
Treses' 8  yield  criterion  Is  also  conipatlble  with  the  Second  Lav  of  Ther¬ 
modynamics  . 

3-5  A  Special  Case  of  the  General  Flew  Rule 

We  new  proceed  to  exhibit  a  further  Justification  of  the  form 
(3.2.8)  of  the  flow  rule.  Specifically^  in  this  section  It  will  be 
shewn  that  the  Von  Mises  flew  rule  formulated  for  perfectly  plastic 
media  Is  a  special  case  of  the  general  flew  rule. 

By  (3.1.10),  the  general  flow  rule  ceua  be  written  as 


df 

“5t“ 


’IJ 


il 


Sh  / 

(  57 


df 


5h  :  V 


mn 


mn 


(5.5.1) 


With  the  yield  function  of  Von  Mises  represented  by 


—  V  V 

2  mn  mn 


(5.5.2) 


where  K  Is  the  constant  yield  stress  In  simple  shear,  the  flow  rule 
now  becoeies 


ij 


lil 


(▼  V  ) 

'  mn  mn 


w 


(5.5.5) 


•  1  /o 

By  definition  (3-1-2),  t  A*  (■jj.  *kl^  '  *  re-«rr»ng 

ing>  {5*5*3)  now  appears  as 


(v  V  ) 
'  an  an 


ij 


<Vki 


1/2 
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(5.5.M 


n1/2 

Prom  (3.5 *2)  follows  that,  during  plastic  flow,  (v^ 
equal  to  therefore 


( 


▼  y__ 

an  nn 

V  «ia 


1/2 

) 


K 


®kl 


(3.5.5) 


where 


§■  ®kl  ®kl 


Combining  (3-5.5)  and  (3-5.^),  we  arrive  at 


Y 


iJ 


K 

/t 


ij 


(3.5.6) 


,  The  e3q>res8ion  given  by  (5-5.6)  is  the  exact  form  of  the  Mlses 
flow  rule  for  perfectly  plastic  media  [3-13].  T^i*  implies  that  the 
genered  flow  rule  as  i^piesented  by  (5-2.8),  under  certain  assuapti'on 
can  he  reduced  to  some  known  plastic  stress  strain  relation,  e.g.  the 


Mlses  flow  rule. 
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PART  n.  AmiCATXOH 

C3Mqpter  4.  Ih«  Hutle-PlABtle  ResjKmse  of  a  Half  Sipaea  to  a  Qalfoxaly 
Apiplled  Heat  Poise  at  its  Boundary 
4-1  mtrodnetlon 

Teaq^erature  variations  of  large  auplitode  and  short  duration  are 
encountered  In  various  fields  of  engineering  applications.  For.exaia- 
ple^  the  elements  of  a  nuclear  reactor  may  he  siibjeeted  to  sudden  heat¬ 
ing  and  cooling  due  to  changes  in  the  rate  of  fission.  Similarly, 
spacecraft  Is  subjected  to  heat  pulses  associated  vlth  the  ignition  of 
rockets,  and  vlth  re-entry  Into  the  atmosphere  of  the  earth  [4.1], 
[4.2],  [4.3]*  It  Is  naturally  Inportant  to  be  able  to  predict  the  tran¬ 
sient  and  residual  stresses  and  defoznatlons  caused  by  such  heat  pulses 
and  to  Investigate  the  physical  damage  Incurred. 

Transient  thermal  stress  and  defoxnatlon  In  the  elastlc^lastlc 
range  have  been  recently  studied  by  a  number  of  Investigators.  In  par¬ 
ticular,  various  problems  for  plates  have  been  treated  by  Weiner  [4.4], 
Tuksel  [4.5]^  Landau  and  associates  [4.6],  and  most  recently  by  Mendel- 
son  and  Bpero  [4.7]  •  Except  for  the  vork  of  Mendelson,  the  Incremental 
theory  has  been  used  In  every  case,  and  the  medium  considered  vas‘ assum¬ 
ed  to  be  elastic,  perfectly  plastic  obeying  the  Von  Mlses  yield  condi¬ 
tion  and  possessing  a  constant  yield  stress.  Landau  [4.8],  however.  In 
one  case  did  take  the  dependence  of  the  yield  stress  on  temperature  in¬ 
to  eonsidermtlon.  Mendelson  employs  the  deformation  theory  vlthoot  re¬ 
sorting  to  any  flow  role  of  yield  condition,  but  takes  strain  hardening 
and  dependence  of  the  physical  properties  of  the  medium  rpon  tempexmr- 
ture  into  account.  Similar  investigations  pertaining  to  elastic,  per- 
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fectly  plastic  cylinders  were  considered  by  Weiner  and  RoddleBtOQ  [4.9]> 
Landau  and  Zirlcky  [4.10].  In  the  latter  Inyestlgatlai  the  yield  condlr 
tlon  of  Von  Mlses  Is  adopted,  and  the  yield  stress  Is  assumed  to  d^end 
on  the  tenqperature.  A  similar  problem  pertaining  to  elastic,  strain 
hardening  spheres  vlth  a  consteuit  yield  stress  vas  treated  by  Huang 
[4.11];  because  of  the  con^lexlty  of  the  numerical  calculatlcsis  In  the 
work  of  Huang  and  Landau,  hl£^  speed  cooputers  vere  used. 

In  order  to  Illustrate  the  use  of  the  concepts  developed  In  the 
first  part  of  this  study,  ve  shcdJ.  present  a  detailed  analysis  of  the 
transient  and  residual  stresses  and  defomatlons  Induced  In  a  Vm-T-P 
space  by  a  heat  pulse  uniformly  distributed  over  Its  boundary.  The 
half  i  i)ace  Is  assumed  to  be  constrained  against  lateral  motion,  euid  the 
medlTsn  Is  assumed  to  be  elastic,  homogeneous  and  Isotropic,  and  linear¬ 
ly  strain  hardening  in  the  plastic  range.  Moreover,  the  yield  stress 
is  permitted  to  vary  linearly  with  temperature,  vhereu  other  material 
properties  ai*e  assumed  constemt. 

We  shall  give  a  cooiplete  stt^  of  the  growth  Md  decay  of  various 
regions  of  loading  and  unloading.  The  most  serious  limitation  of  the 
solution  to  this  problem  Is  that  It  lacks  a  characteristic  dimension  of 
length.  Therefore,  we  shall  briefly  Investigate  a  plate  of  finite 
thickness  In  Chapter  6. 

4-2  The  Temperature  Problem 

As  a  consequence  of  the  mathematical  cosplexlty  Inherent  In  the 
thermal  stress  problems,  it  Is  necesssury  to  select  teiqperature  solu¬ 
tions  that  are  as  slsple  as  possible.  In  peurtlcular,  the  mathematical 


problew  aore  greaitly  slnpllfled  if  the  tem>ex«ture  eolation  is  In  a 
eloeed  fozn.  One  class  of  such  closed  foxn  solatitnui  corresponds  to 
tesqperature  fields  Induced  by  Instantaneous  heat  sources  (or  sinks), 
doublets,  and  the  like.  Since  these  solutions  are  singular  at  the 
source  point,  the  latter  must  be  located  outside  the  material  body. 


Fig.  4-1  Heat  Sources  and  the  Infinite  Slab 

In  the  present  analysis,  a  row  of  beat  sources  is  placed  at  an 
arbitrary  distance  a  from  the  boundary  of  the  body  (Fig.  4-l).  Let 
X  be  the  coordinate  meaexured  from  the  source  plane,  and  t  the  time. 
The  Increewe  of  teBg>erature  T  over  seme  reference  teo^erature  T^  Is 
then  given  by  [4.12] 

-  x^/lwct 

T  -  C  -2^  (4.2.1) 

vhere  C  is  a  constant,  and  K  is  the  dlffiisivlty  of  the  material 
(assumed  to  be  constant).  It  nay  be  readily  verified  that  the  tempera¬ 
ture  field  (4.2.1)  satlsifes  the  one-dlmecsloneil  heat  ccmductlon  equa- 
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tiCQ 


i? 


1  a® 

K  5t 


ik.2.2) 


and  Is  finite  eveiTvhere  In  the  material  (l.e.  a). 

The  ataximum  ten^erature  T  occurs  at  the  boundary  x  -  a.  and, 

m 

as  may  be  fotmd  from  (4.2.1),  at  the  Instant  t  given  by 

n 


t  -  2- 
m  ~  2k 


(4.2.3) 


Similarly  it  is  found  that  the  maximum  temperature  at  x  is  reached  at 

the  Instant  t  given  by 
mx  “  •' 


mx 


2 

X 

2k 


(4.2.4) 


The  arbitrary  constant  C  in  (4.2.1)  may  now  be  eliminated  in 

terms  of  a  given  maximum  temperature  T  ,  thus  obtaining 

m 


~  -  X^/kKt 


m 


(2  K  t) 


W 


(4.2.5) 


Substituting  (4.2.4)  into  (4.2.5),  the  maximum  tesqperature  T 


mx 


attained  at  x  becomes 


T  -  T  - 
mx  m  X 


(4.2.6) 
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FIG  4-2  TEMPERATURE  VARIATION  WITH  TIME  AT  DIFFERENT  POSITIONS  :(4.25) 


39 


FIG.  4-3  TEMPERATURE  DISTRIBUTION  AT  DIFFERENT  INSTANTS  OF  TIME 
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She  vazlatloB  of  the  taqperature  T  with  tlae  t  for  'various  -values 
of  position  X  is  shewn  In  Tig.  4-2,  whereas  Fig.  4-3  illustrates  the 
-variation  of  T  vl-th  x  for  -various  -values  of  t.  nie  results  ha-ve 
been  made  dinenslatiless  by  using  the  -variables 

Va,  t/t^ 

for  ten^rature,  position  and  time  respecti-vely. 

The  solution  (4.2.1)  corresponds  to  a  teoperature  field  Induced 
by  suddenly  releasing  a  finite  amotmt  of  heat  at  x  s  0.  It  may  be 
seen  from  Fig.  4-2  that  -this  solution  is  a  reasonable  approjdmatlon  of 
a  heat  shock,  l.e.,  of  a  sudden  rise  in  the  boundary  teoperature  follow¬ 
ed  by  a  gradual  decay.  It  has  been  asstnaed,  of  course,  that  -the  heat 
conduction  problem  is  not  influenced  by  -the  mechanical  defozmation  of 
-the  material. 

4-3  The  Elastic  Regime 

During  the  Initial  s-tages  of  -the  heat  shock,  the  response  of  -the 
material  will  be  elastic.  Let  €^y  u^  denote  stresses, 

strains  and  displacements  respecti-vely.  The  elastic  thermal  stress 
problem  is  -then  characterized  by  the  stress  eqiiatlon  of  equilibrium 


0 


The  stress  strain  relations 


(4.3.1) 


Bi  a  T  -  (1+1/,)  5,^ 


(4.3.2) 


kl 


and  the  strain  dlsplaeenent  relations 


5  (“l.J  ^  "j.!* 


(l^.3.5) 


yhere  a,  aure  Young's  aodulus,  coefficient  of  thexnal  expan¬ 

sion  and  Poisson's  ratio,  respectively,  and  use  has.  been  made  of  the 
eonventianal  index  notaticm. 

In  addition  to  quiescent  initial  conditions,  it  is  assumed  that 
at  the  boundary 

V 

T  ■  0  for  X  =  a  (4.3A) 

In  view  of  the  simplicity  of  the  problem,  the  following  abbreviated 
notation  is  Introduced, 

U  -  U  (x,t)  ,  =  U„  -  0  (4.3.5) 

X  jr  « 

€  =  €  (x,t)  ,  e  =  e  ■  0  (4.3.6) 

XX  '  '  '  *  yy  zz  \  '  t 

T  ■  T  ■  T  (x,t),  T _  ■  0  (4.3.7) 

yyzz'^'^  XX  \  ^  m 

The  relations  (4.3.5)  es^ress  the  condition  that  the  hadf-space  is  con¬ 
strained  aigainst  latere^,  motion.  It  follows,  then,  that  adl  sheaxing 
stresses  and  shearing  stralps  vanish.  Finally,  the  second  of  (4.3.7) 
follows  froB  (4.3.1)  and  (4.3.4). 

Taking  (4.3.5),  (4*3.6),  (4.3.7)  into  alccount,  the  relation 
(4.3.2)  finally  aqppeanrs  as 


k2 


Ej^€-E^aT«-2l/j^T  , 

-a  T  -  -v  (1  -v^)  / 


whence  follows 


e  s  e 


XX 


1  + 
1  - 


a  T 


(4.3.8) 


a 


®l 


"  -  V  -  -  -  1  > 


(4.3.9) 


Here  T  is  given  by  (4.2.5),  €md  we  note  that 


T  - 


(ll.S.lO) 


where 


\ 


2(1  +  \) 


(4.3.11) 


Is  the  elastic  modulus  of  shear.  Introducing  the  maximum  shearing 
stress  q. 


a  E, 


2(1  -  74) 


(4.3.12) 


and  the  maximum  shearing  strain  p 


P  -  2  - 


€ 


2  (1  -  \) 


a  T 


(4.3.13) 


1^3 


th«  relation  (4. 3  >10)  ney  nor  be  written  at 

<1  -  2  Ug  p  (4.3.14) 

Finally,  Integrating  (4.3. 8)^  the  dljqplacaBent  U  It  obtained  at 

1  +  ^1  *  V2 

U  - -  a  a  T  (e  |)  I  if  ( — -  ij  (4.3.15) 

1-Vj^  ®  a/TT 

where  , 

8  a 
a  r  ■ 

“  7^/a  j  ® 

Is  the  error  ftmctlon,  and  the  arbitrary  function  of  time  hat  been 
choten  to  that  the  regularity  condition 

U  -  0  for  X  -♦  -o  (4.3.16) 

It  satltfled. 

Ve  also  make  thd  relevant  observation  that  the  solution  given  by 
(4.3.8),  (4.3.9)^  (4.3.15)  holds  even  If  the  elastic  region  begins  at 
tone  elastic  plastic  Interface  characterized  by  x  ■  p^. 

4-4  Pie  Elaatlc -Plastic  Regime 
4-41  !nie  Onset  of  Yielding 

It  Is  found  from  the  elastic  solution  presented  In  the  preceding 
section  that,  at  the  teoQ>erature  Increases,  the  shearing  stress  also  In- 

creates.  Let  the  Initial  yield  stress  In  sheso*  be  denoted  by  y^,  and 
let 


yj  -  -  HP 


kk 


vhere  B  Is  assumed  to  loe  a  knoim  constant,  and  Is  th^  value  of 
y^  for  T  ■  0  (l.e.  at  the  reference  teiqperature  T^).  The  Instant 
t^  at  which  yielding  begins  at  the  boundary  x  «  a  can  be  readily  de¬ 
termined  by  letting 

4  (a,t)  -  yj  (a,t)  ,  (a) 

if  we  adopt  the  yield  condition  of  Tresca.  From  (4.3«12),  (4>4.l)  and 
(4.2.5)  we  then  obtaija 

loe*!*  rrt^  ■ 

(4.4.2) 

The  equation  (4.4.2)  can  also  be  used  to  find  the  instant  t^  at 
which  yielding  occurs  at  the  plane  characterized  by  x,  if  a  and  t^ 
are  replaced  there  by  x  and  t^,  respectively. 

From  (4.5.12),  (4.4.1)  and  (a)  follows  that  the  boundary  teiapera- 
ture  at  the  Instant  t^  is  given  by 


yp 

B  +  a/2(l  -  V^) 


(4.4.3) 
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Fig.  4-4  Regions  of  Loading  During  the  KLastle-Flaetlo  Reglae 

For  t>' an<  elastic  pleistlc  Interface  whose  location  Is  denot¬ 
ed  hy  progresses  into  the  material.  In  the  vake  of  this  inters 

face,  the  material  uidergoes  plastic  loadixig  until  t  becomes  eqoeil  to 
t^;  then,  corresponding  to  the  decrease  of  the  boundary  temperature, 
an  elastic  unloeding  front  emanates  from  the  boundary. 

The  present  section  concems  the  time  Interval  t^  ^  t  ^  t^,  hence. 
In  addition  to  the  solixtlon  of  the  elastic  region,  it  Is  necessary  to 
derive  a  solution  for  the  plastic  region. 

4-42  Ihe  Elastic  Region 

The  solution  of  the  elastic  region  Is  the  same  as  that  given  in 
section  4-3  except  that  the  elastic  region  now  lies  at  x>  p^,  where 
p^  Is  yet  to  be  determined. 


In  order  to  detemlne  the  position  of  the  elastic  plastic  inter¬ 


face  ve  nse  (4.2. 5)>  (4.4.1),  (4.3.12)  and  (4.3.9)*  In  partleu- 
lar,  f]^ 


4  (pjL>  ^^1^*^ 


(4.4.4) 


follows 


1  fli 

2  1-Vj_ 


T  (Pl,t) 


-  B  T  (p^,t) 


or 


(4.4.5) 

It  is  Interesting  to  note  that,  m  a  consequence  of  (4.4.4),  the  ten- 
petature  at  the  Interface  remains  constant  and  equal  to  (the  bound¬ 
ary  tenq>erature  corresponding  to  incipient  yielding),  thus 


T(pj^  t)  - 


B  +  a/2(l-1/^) 


(4.4.6) 


where  is  the  tenq>erature  corresponding  to  Incipient  yielding  at 

the  plane  characterized  by  x. 

It  follows  then,  that  the  progress  of  the  interface  my  be  studied 
by  drawing  on  Fig.  4-2  a  horizontal  line  at  the  hel^t  of  Aod 

investigating  its  Intercepts  with  tenperattire  cuirres. 
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Since  la  the  maylmm  tegqperature  for  yielding,  the  aaxlna 
value  of  can  he  detemlned  eaally.  Bunely,  since  the  waximua 
teBqperature  T  reached  at  x  Is  given  hy 


(4.2. 6)R 


It  f  oUcsrs  that 


or  hy  (4.4.6) 


(4.4.7) 


Therefore,  It  Is  known  that  the  region  a  £  x  <  viU  xindergo  sane 
plastic  defonaatlon.  The  exact  detezninatlon  of  the  residual  state, 
however,  can  only  he  found  frcm  a  detailed  analysis. 


The  Plastic  Region  of  Loeidl 


As  the  houndary  temperatme  rises  heyond  T^,  the  pleustlc  region 
extends  Inward  from  the  boundary.  In  order  to  derive  the  appropriate 


solution  for  the  plastic  region,  the  plastic  stress  strain  relations 
must  he  estahllshed  first;  this  step  must  he  preceded  hy  the  formula¬ 


tion  of  yield  functions,  hcwever. 

4-45-1  The  Yield  Functions 

In  the  elastic  solution  the  maximum  shearing  stress  g  and  maxi¬ 


mum  shearing  strain  p  la  related  hy 
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4  -  2  Hj.  P  (4.3.14)r 

since  the  medium  Is  assvimed  to  be  linearly  strain  hardening,  the 
relation  between  the  maximum  shearing  stress  and  maximum  shearing  strain 
Is  of  the  foxm  Illustrated  graphically  In  Flg«  4-3  • 


Fig.  4-3  Relation  Between  Maximum  Shearing  Stress  and  Strain 
In  a  Linearly  Strain  Hardening  Medium 

We  write 

yc  =  ^  y  " 

where  y  Is  the  c\aTent  yield  stress,  y_  Is  the  initial  yield 
stress  which  Is  equal  to  y^  -  BT,  emd  ^  Is  the  Increase  In  yield 
stress  due  to  strain  hardening.  Similarly,  we  let 

dp  »  p"  +  A  p'  ,  (b) 

where  Is  the  Increase  in  total  strain  corresponding  to  the  In¬ 
crease  in  yield  stress,  p"  Is  the  plastic  component  of  the  maxliinB 


■hearing  strain,  euid  Ap*  Is  the  corresponding  Increase  In  the  elastic 
conponent  of  the  maximum  shearing  strain. 

Referring  to  Fig.  4-5,  we  note  that 


Ay  «  2  ^Jlg  ,  (c) 

and 

Ay  =  2Wp(p"  +  Ap‘)  (d) 


We  eliminate  Zip'  from  (d)  hy  combining  (c)  €uid  (d). 


Ay 


1 

1-s 


2  p  p" 

P 


where 


s 


(e) 


From  (a)  and  (e),  the  current  yield  stress  Is  then  obtained  In  the 

form 

Since,  by  (4.5.15),  the  maximm  elastic  shearing  strain  p'  Is  de¬ 
fined  as  £^05/2,  we  have  that 


Zp'  -  A  £30/2 


(g) 


Similarly,  we  write 
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vhere  ia  the  plastic  strain  coeq;>oaent  in  the  z  direotioa.  It 
rosains  to  esqpress  p"  in  terns  of  the  strain  hardening  paraaeter  ( 
defined  hy  (^.1.2). 

leaking  the  plastic  Inccopresslblllty,  and  the  symmetry  about  the 
X  axis  into  account «  ve  rec^lly  find 


S  B  S  B 

yy  az 


’  2  ®xx  ^ 


Moreover, 


S  bS  bS  bO  , 
xy  yz  xz  ’ 


and,  by  (3.1.2), 


r  2  1  2  V2 

S  ■  j  -  I  'xx!  '  “ 

°  1/2 

-  <l>  'xx 


(1) 


MsMng  use  of  (h)  and  (l),  ve  rewrite  (f)  as 


yc 


BT  + 


(J) 


If  the  medium  has  not  undergone  previous  plutic  deformation,  the  para¬ 
meter  I  is  equal  to  zero,  (j)  reduces  to  (It-. 4.1). 

4-45-2  The  Stress  emd  Strain  Rates  Corresponding  to  Positive  Testpera- 
ture  Rate 

The  yield  criterion  of  Tresca  will  be  used  here  not  only  to  pro¬ 
vide  a  means  for  locating  the  elastic -plsistlc  Interface,  but  atlso  to 
serve  as  yield  function  from  which  plastic  stress  strain  relations  are 
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darlTed.  Ve  observe  froei  the  elastic  solution  that  the  stresses  t _ 

yy 

and  T  .  are  nesatlve  corresponding  to  rising  tsogperature,  and  posl> 
zs 

tlve  for  falling  teiqperattnre.  The  corresponding  maxinam  shearing 
stresses  are,  therefore,  given  by  -r^y/g,  1“  first  case, 

and  by  ‘‘’yy/S*  ^  second  case. 

Ve  now  ooinblne  (3.3*3),  (3*3.^)  vlth  (j),  euid  obtain  a  non-lso- 
thezml  yield  condition  in  the  form 


-  H  -  -  (y^  -  BT  +  ~  )  (1|..4.8) 


1-6  ^p 


r(2)  .  ,<2)  .  ,  .  js^  -(y^-ar.  ^ 


)  ih.k.9) 


Where 


,(i) 


V  -  V 

JE _ L 

2 


Ijol 

2 


.(2) 


V  -  V 
X  z 


zz 

2 


(4.4.10) 


H  -  r,  -  (y„  -  HT . 


In  particular,  since 


T  »  T  «  T 

yy  zz 


ve  have  that 


f(l)  ,  f(2)  ,  f  .  .  t/2  ,  -  F^®^  «  F  (4.4,11) 
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Fig.  k-6  Yield  Surface  and  the  Path  of  Loading  for  the  Present  Prohlem 


Prom  (4.4.11)  follows  that  the  state  of  stress  In  the  plastic  reg¬ 
ion  corresponds  to  the  yield  comer  which  is  represented  by  the  point  A 
In  Fig.  4-6.  The  plastic  stress  strain  relations  may  thus  be  obtained 
from  (5.3.1),  (3.3.2),  and  (3.2.8): 


(7) 


ij 


7-1 


•iJ 


(7) 


’iJ 


_ Ij _ 


kl  "^kl 


(4.4.13) 
+  ap/Sr  T) 


(y) 

where  f  are  given  by  (4.4.8),  (4.4.9)  uid  (4.4.10).  The  flow  rule 
(4.4.13)  Is  ccnplemented  by  Hooke's  Law,  (4.3.2),  and  the  condition  that 
the  total  strains  In  the  y  and  z  directions  vanishes.  These  condi¬ 
tions  furnish  five  equations  for  the  five  unknowns  £  .  s _ ,  £ 

XX  XX  yy 
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■ _  and  T  in  the  plastic  region.  3peclfleally,  we  find* 

77  77 


s  « 

XX 


~L 

1-.  “p  ^  3  ' 


( -I 


(») 


a  m  B  m 

yy  zz 


1 

2 


Z 

l-«  '^P  ^  3  ^ 


(  -I  (^) 


£  »  2  +  O  T 


XX 


tt 


^yy  V  ^  ^ 


(c) 

(d) 


8  “  “  £ 
yy  ‘^yy 


(e) 


Since  the  stress  t  =  t  _  and  the  total  strain  e _  are  of  prl- 

yy  zz  joc 

nary  Inportance,  ve  solve  for  then  front  (h.4.1h)p  and  obtain 


yy  zz 


M  T 


XX 


«  +  £  ^ 
XX  XX 


N  T 


(4.4.16) 


fhe  derivation  of  (a)  and  (b)  In  (4.4.14)  Is  given  In  Ai^ndlx  B. 


tfhere 


1/2 

2  BEj^  (l-s)  -  2  g  Up  W) 

(1-8)  +  2  (1-1/j^)  Up  (4/5)^^^ 


(4.4.17) 


1/2 

H.  fa.  Bjl^ilil) - j 


J 


Ve  ai^haslze  that  equations  (4.4.15)  and  (4.4.16)  are  applicable 
only  to  the  plastic  regions  where  plastic  flew  Is  Induced  by  rising  tem¬ 
perature.  Moreover,  with  reference  to  (4.4.15)  and  (4.4.1?),  we  note 
that  a  positive  temperature  rate  nay  give  rise  to  a  positive  or  nega¬ 
tive  stress  rate  depending  on  whether  the  effect  of  tenq>erat\ure  upon 
yield  stress  domlmtes  or  the  effect  of  strain  hardening  upon  yield 
stress  dominates.  In  theflrst  case,  the  stress  t  In  the  plaatlc 
region  will  always  be  less  than  the  valtie  It  aBB\]ned  at  Inltlskl  yield¬ 
ing. 

The  stress  and  strain  distribution  In  the  plMtlc  region  may  be 
obtained  by  Integration  of  (4.4.15)  and  (4.4.16). 

4-43-3  The  Stress  and  Strain  Bates  Corresponding  to  Negative  Tempera¬ 
ture  Rate 

If  the  boundai^  of  the  hali*  space  Is  cooled,  the  stress  in  the  y 
and  z  directions  becomes  tensile  emd  the  maximum  sheafing  stress  be¬ 
comes  Tyy/2  ■  ‘fjjg/S  in  the  place  of  -  t^^/2  •  -  f^^/2  found  for  the 
case  of  heating.  The  yield  functions  corresponding  to  (4.4.8)  amd 
(4.4.9)  new  assume  the  forms 


55 


,(!)•  .  ,(!)• 


HP  + 


P  /T 


p(2)'  .f(2)' 


-  HP  + 


(k,k.ia) 


(4.4.19) 


vhere 


f(2)' 


H 


y  -  HP  + 


1 

1-8 


2  K 


/r 


A^ln,  since 


T  =  T 
yy  zz 


ve  have  tbat 


,(i)'  .  ,(a)' 


f'  ,  .  p 


Proceeding  as  'before.  It  Is  found  tbat  the  stress  rate  and  strain 
rate  are  related  to  the  ten^erature  rate  by 


(4.4.20) 


€  .  H'T 

XX 


(4.4.21) 
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vhere 


M' 


1/2 

-2  B  Ej^(1-8)  -  2Ej^  a  u  (4/5) 

1/i  " 

(l-a)  +  2  tip  (1-14)  (4/^) 


(4.4.22) 


1/2  1/2  2.V 


P 


Relations  (4.4.20),  (4.4.21),  (4.4.22)  nay  only  be  used  to  describe 
the  stress  and  strain  rate^  in  any  plastic  regions  vhere  plastic  flew  is 
induced  by  falling  teB5)erat\ire . 

If  the  j'leld  stress  is  independent  of  temperature,  the  response  of 
the  medltmi  should  be  the  same  regeordless  of  whether  the  external  excita¬ 
tion  is  a  heating  or  a  cooling  pulse.  This,  in  fact,  is  the  case  if 
the  quantity  B  in  (4.4.15)  and  (4.4.20)  is  set  equal  to  zero,  in  which 
case  these  equations  become  identical;  the  same  is  true  for  (4.4.16)  and 
(4.4.21).  Purthennore,  if  in  addition  to  setting  B  equal  to  zero  in 
(4.4.15),  (4.4.20),  (4.4.16)  and  (4.4.21),  the  plastic  shear  modul^ls 
is  set  equal  to  the  elastic  shear  modul\is  these  equations  re¬ 

duce  to  their  elastic  counterparts  represented  by  (4.3.8)  and  (4.3.9) • 
4-43-4  The  Criteria  of  Loading  and  Ihaloading 

As  was  noted  previously,  the  medium  will  undergo  scsie  plastic  de- 
fomatlon  in  the  region  a  $  x  ^  R^  if  temperature  exceeds  the  value 
T^  given  by  (4.4.6).  For  a  point  in  this  region,  two  instants  axe  of 
prime  importaurce:  the  instant  at  which  plastic  flow  begins,  and  the  in¬ 
stant  at  which  plastic  flow  ceeuses.  The  first  instauit  can  be  dbtalned 
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frcm  (4.4.2),  vhereas  the  second  Instant  renains  to  he  detexnlned. 

Hhe  loetdlng  and  mloadlng  criteria  discussed  In  Chapter  3  nay  helj> 
to  shed  some  ll^t  on  the  duration  of  plsistlc  flew.  From  the  relations 
(5.5.2),  (4.4.8),  (4.4.9),  (4.4.11),  (4.4.12)  follows  that,  for  the 
present  case,  these  criteria  may  he  esqpressed  hyt 


-i  Bi>0 

-i  m<0 


for  loading 


for  unldading 


(4.4.25) 


Substituting  (4.4.15)  into  (4.4.25),  the  cirteria  become 
(B  -  |m)  T>0  for  loading 

(4.4.24) 

T  * 

(B  -  |M)  T  <0  for  unloading 


where  M  is  given  by  the  first  of  (4.4.17).  The  quantity  inside  the 
parentheses  of  (4.4.24)  is  a  constemt,  emd  is  found  to  be  positive 
definite  upon  closer  examination.  Therefore,  the  criteria  for  loading 
euid  unloeiding  may  be  stressed  as 


T  >  0 

T  <  0 


for  loetdlng 
for  vinloadlng 


(4.4.25) 


t  The  derivation  is  given  in  Appendix  B. 
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It  Is  thus  seen  from  (4.4.2^)  that,  at  any  location  in  the  region 

represented  hy  a  <  x  ^  plastic  flow  begins  at  the  instant  t^ 

at  which  the  teng>erature  T  is  at  emd  temlnates  at  the  Instant 

t  „  for  which  the  temperature  is  at  its  maximum  T  . 
mx  mx 

4-43-5  Solution  for  the  Plastic  Region 

The  stress  T  and  the  total  strain  e  in  the  plastic  region 
yy  XX 

may  be  obtained  by  integrating  (4.4.15)  and  (4.4.16)  between  the  limits 
tj^  and  t  where  t^  is  the  instsuit  of  incipient  yielding  at  x. 

Rjr  (4.3.8),  (4.3.12)  and  (4.4.6),  we  obtain  after  appropriate  integra¬ 
tions 


M  T(x,t)  -  [2  y„  + 


B  +  ^a/2(i-i/^; 


(M  -  2B)] 


(4.4.26) 


«xxf*»*^  =  W  T(Xit) 


r/-  ^*'^1  ... 

[(N  -  1  ®)  B  +  a/2(l  -  T/^)  ^ 


(4.4.27) 


where  M  and  N  are  given  by  (4.4.17). 

The  displacement  U^,  denoted  by  U,  is  obtained  in  the  plastic 
region  by  integrating  (4.4.27)  with  respect  to  x  between  the  limits 
and  X 


1+74 


u(x,t)  ^  (H  -  (Pl  "  *) 


+  g  (P^^t) 


9 


(4.4.28) 
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vhere  tCs)  is  the  error  function  defined  by  (1|-.3»15)*  The  function 
osy  be  determined  by  ccsislderlng  the  condition  of  continuity 
of  displacements  at  the  elastic  plastic  Interface  x  ■  namely 


U 


eleuBtlc 


U. 


plastic 


at 


(4.h.29) 


Prom  (4.4.28),  (4.4.29)  and  (4.3.15)  then  follows  the  result 

^  \  (  I  tt(-4zr)  -  1] 

1  2^/7Ct 

The  congplete  solution  for  the  time  interval  t,  ^  t  ^  t  is  given 

1  m 

by  (4.3.8),  (4.3.9)>  (4.3.15)  for  the  elastic  region  and  by  (4.3.26), 
(4.4.27),  (4.4.28),  (4.4.29)  for  the  plastic  region. 

4-5  The  Bastlc -Plastic-Elastic  Regime 
4-51  The  Onset  of  Elastic  Dhloadlng 

This  section  Is  concerned  with  the  response  of  the  material  at 
time  \  where  t^  is  the  instant  at  which  the  temperature  at  the 
boundary  is  at  its  maximum.  We  recall  here  the  criteria  (4.4.25)  for 
loading  and  unloading,  and  observe  that  any  point  in  the  region 
a  rsx  unloading  is  imminent  when  T  becomes  zero.  In  other  words 

unloeidlng  Inpends  as  soon  as  the  temperature  at  any  location  in  the  re¬ 
gion  a  <  X  <  has  attained  its  peak  value . 

Prom  inspectiwi  of  (4.2.3)  and  (4.2.4)  it  is  found  that  the  tem¬ 
perature  first  attains  its  maximum  at  the  boundary  at  the  instant  t 

ffl 

g^-ven  by 


(4.2.3)R 
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At  this  very  Instant  the  teniteratitre  at  the  boundary  begins  to  decrease, 
and  an  elastic  unloading  front  starts  to  progress  Into  the  medium.  The 
unloading  front  eventually  overtakes  the  front  of  plastic  loading  at 
the  Instant 


(4.2.4)R 


and  at  the  position  x  =  of  the  half  space. 


1 

1 

p 

x=R- 

Region 

^Region  of 

A  Elastic 
y  tJhloedlng 

Vy 

1 

Region  of 
Plastic 
Loading 

"1 

Elastic 

_ X»D-  _ ^ 

First  Elastic 
Plastic 
Interface 

1 

Second  Elastic 
Plastic 
Interface 

Fig.  4-7  Regions  of  Loading  and  Uhloedlng  In 
the  Elastic -Plastic -Elastic  Regime. 


4-52  The  Elastic  Region 

For  time  t  >  t  ,  the  elastic  region  still  lies  beyond  the  elastic 

SI 

plastic  interface  the  position  of  which  Is  denoted  by  and  may  be 

found  from  (4.4.5).  The  solution  for  this  region  is  identicsJL  to  that 
derived  In  sectlcsi  4-?. 


6l 


4-53  The  PlMtlc  RegKm  of  Loading 

The  plastic  region  now  lies  between  the  two  Interfaces  and 
where  pg  denotes  the  position  of  the  elastic  unloading  front. 

In  this  region  the  equatlc»is  (4.4.2^),  (4.4.27)  *uid  (4.4.28),  derived 
In  the  last  section,  still  represent  a  valid  solution. 

Let  Tis  recall  that  the  unloeuiilng  front  Is  located  at  a  point  where 
the  ten^emture  Is  at  Its  saxlmum.  For  a  given  position  x,  the  In¬ 
stant  t  at  which  the  ten^erature  attains  Its  maxlmuni  Is  found  from 
ax 

(4.2.4)  to  be 


t  =  x^/2K  (4.2. 4)R 

nx 

Therefore,  the  position  x  =  Pg  of  the  unloading  front  Is  given  by 

p/  =  ac  t  ,  (4.5.1) 

5  t  >  t  , 

X  m 

where 

=  Rj^/2k  ,  =  a^/2>c 

The  position  x  =  P^  of  the  plastic  loeidlng  front  can  be  determined  as 
befbre  from  (4.4.5)* 

5-54  The  Elastic  Region  of  Unloading 

For  any  Instant  t  between  t  _  and  t  ,  three  regions  exist 

M; 

simultaneously  In  the  half -space,  namely,  the  elastic  reglcm  of  unloeid- 
Ing,  the  plastic  region  of  loeidlng,  emd  the  elastic  region.  Hcwever, 
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tor  i  >  only  the  elastic  region  of  unloading  and  the  elastic 

region  renaln.  R^laclng  x  hy  In  (h.2.4),  the  Instant  t^ 
may  he  obtained  by  substituting  the  e^^resslon  for  R^  given  by  (4-.  4.7) 
In  the  resulting  equation,  thus. 


2(l-Vl) 


(4.5.2) 


4-55  The  Transient  Solution  In  the  ELeustlc  Region  of  ttiloeMilng 

The  stress  and  strain  rates  In  the  elastic  region  of  unloading  nay 
be  obtained  from  (4.5.8)  and  (4.5.9) 


T 


(x#t) 


E^a 

1^ 


i  (x,t) 


(4.5.5) 


1+14 

€  (x,t)  =  ®  ^  (x>t) 


(4.5.4) 


Ve  shall  take  as  the  limits  of  Integration  t  Euid  t  ,  vhere 

nx 

t  Is  greater  than  t^.  Ve  also  recall  that,  by  (4.2.6),  the  maximum 
teDq>erature  T(x,  t^)  is  equBd.  to  Then  tislng  for  t(x,  t^), 

€(x,  t^)  the  results  (4.4.26),  (4,4.27),  integration  of  (4.5.3), 
(4.5.4)  yields 


T(x,t) 


*  [  (  ^^  *  *)  I 

for  t  >  t 
mx 


(4.5.5) 
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€(x,t) 


1+Z4  1+24  .  y_ 

1^  “  7  "  B  +  E^a/2(l-2^p^ 

for  t>t  .  (^•5»6) 

mx  \  ✓  / 


In  order  to  obtain  the  displacement  U,  we  Integrate  (4.5.6)  between 
the  limits  and  x  and  obtain 


l+K  ^  1/2  /_x _ V  __  /  ^1 

U(x,t)  =  a  a  j  e)  ^^^2^ 


1+K  V 

"l-U^  ^  ^  "  B  +  a/2(l-J/j^;  ■  ^1 


)] 


+  g  (pi>  t) 


(4.5.7) 


where  is  defined  eis  follows; 


P 


1 


f 


for  t  _>  t  >  t 
mRj^  n 

for  t  >  t  _ 
mR^ 


Fozmally  the  function  s(pj^/  t)  in  (4.5.7)  is  deteimined  from  the  con- 
dltlcxis  of  continuity  of  displacements  across  the  interface  x  = 


^elastic  unloading  “  Aplastic  * 
at  X  .  Pj  for  t  >  t.  , 


(I..5.8) 
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and 


^elastic  vtnloadlng  “  ^elastic  , 


at  X  =  R,  for  t  ^  t  _ 
X  xnx\« 


ThTis,  by  (4.5.7),  (4.5.8),  (4.5.9),  (4.4.a8),  (4.4.30)  and  (4.3.15), 


g(p2>t)  =  (»  -  Cl)  B  +  a/2(l-V^;  (^1  ‘  ^2^ 


■ff  V2  P,  Pp 

N  a  T  (  I  e)  [♦(— ^)  -  ♦(—!:)] 


1+i/,  ^  1/2  p- 


(4.5.10) 


and 


1+U 


1/2 


g(Rl,t)  =  a  a  (  I  e)  [t(— p)  -  1]  (4.5.11) 

X  2 


4-56  steady  State  Solution  In  the  Elastic  Region  of  Unloading 

Expressions  for  residual  stress  t  ,  strain  e  and  permanent 

A  R 

defonnation  Ujj  may  new  be  derived  frem  (4.5.5),  (4.5.6),  (4.5.7)  and 
(4.5.11)  by  letting  the  time  t  in  these  equations  tend  to  Infinity. 
We  readily  find  that 


0®, 


-  [2  y„  + 


o  B  +  ^/Z^l-V^) 


(M  -  2B)]  , 


(4.5.12) 
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i+K  .  y. 

6jj(x)  »  [Tjij  -  -  B  +  a/2{l-p^)^  (^•5.13) 


1+i/.  „  y. 

Up(x)  =  (N  -  3;:5^«)  ^  "  B  + 


1+2/  1/2 


(J*.5.14) 


where,  by  (4.4.7),  the  qioantity  is  given  by 


(4.4,7)R 


We  observe  from  (4.5.12),  (4.5*13)  and  (4.5.14)  that  the  residual 
stress,  strain  snd  deformation  decrease  from  their  maximum  values  at 
the  boundary  to  zero  at  x  =  The  last  term  in  (4.5.14)  represents 
the  uniform  e;Q>anslon  of  the  material,  euid  is  due  to  the  assumed  van¬ 
ishing  of  displacement  at  infinity*. 

With  the  trsuisient  solution  given  by  (4.5*3),  (4*5*6),  (4.5.7)» 
and  the  steedy  state  solution  given  by  (4.5*12),  (4.5*13),  (4.5*14), 
solution  for  the  elastic  region  of  unloading  is  completely  determined. 
However,  it  should  be  noted  here  that  the  relations  Just  mentioned  above 
are  applicable  to  any  position  in  the  region  of  unloewiing  only  if  the 
material  there  is  in  the  elastic  state. 

4-6  Extension  of  the  Solution  to  Higher  Boundary  Temperature 

In  the  preceding  sections  of  this  chapter,  it  was  tacitly  assumed 
that  the  maximum  boundary  tenperature  T^  is  of  such  a  magnitude  that 


*cf.  The  equations  (4.5*15)  and  (4.3*16). 
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no  yielding  occurs  in  the  region  of  unloeiding.  Hovever,  an  inspection 
of  the  equation  (4. ^.12)  for  residual  stresses  shows  that  the  residual 
stresses  Increase  with  Increasing  values  of  the  maximua  boundary  tem¬ 
perature  T  .  It  Is,  therefore,  possible  that  for  stifficiently  large 
m 

T  ,  the  process  of  elastic  unloeuilng  may  be  followed  by  plastic  flow. 

ffl 

In  what  follows  this  possibility  will  be  fully  Investigated. 

4-6l  Yielding  In  the  Region  of  Elastic  TJhloetdlng 

As  was  noted  in  section  4-4,  unloeiding  at  the  point  x  begins  at 

the  Instant  t  when  the  tenroerature  there  has  reached  Its  meuclmum 
mx 

value.  Let  us  new  define  t^^  to  be  the  instant  at  which  plastic  flow 
occurs  during  unloading  at  x.  If  we  account  for  the  temperature  de¬ 
pendence  of  yield  stress,  then,  as  may  be  seen  from  Pig.  4-8,  yielding 

in  the  opposite  sense  will  Impend,  if  the  stress  t  or  t  „  has 

yy  zz 

changed  by  the  amount 

2/t  (x,t  )/+2B[T(x,t  )  -  T(x,t  )] 

'  yy  '  mx  '  mx'  '  '  px" 

Equivalently,  at  the  instant  of  incipient  yielding  in  opposite  sense, 
the  maximum  shearing  stress  q(x,t}  has  changed  by  the  amount 

q(x,t)  =  +  B  [T(x,t^)  -  T(x,tp^)] 

during  the  time  interval  from  the  start  of  unloading  at  t^^  \intil 

the  instant  t  .  The  mathematical  expression  for  the  criterion  of 
px 

yielding  in  the  region  of  elastic  unloading  is  then  as  follows 


FIG.  4-8  PATH  OF  LOADING,  UNLOADING  OF  A  STRAIN  HARDENING  MATERIAL 
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T  (x.t  )  -  T  (x.t  ) 
yy'  *  px'  yy'  '  mx' 


If  equations  (4.2.6),  (4.4.26)  and  (4.3*5)  are  substituted  into 
(4.6.1),  ve  obtain 

E  a  E  a 

(2B  -  T(x,t^^)  «  (23  -  -  2M)  I 

+  2  [2  y^  +  B  +  (4.6.2) 


It  is  of  interest  to  detemine  the  meuclmum  boundary  temperature 

for  which  the  elastic  unloeiding  will  not  be  followed  by  yielding. 

Setting  T(x,t  )  in  (4.6.2)  equal  to  zero,  we  find  the  result 
px 


1  ^0 

^  “  [M  +  QE^/2'a-1K  )  -  bJ  ^^^o  B  +  a/2(l-Z<  ) 


(4.6.5) 


If  the  maximum  boundary  tenperature  is  greater  them.  plastic 

flow  will  occur  again  during  unloading.  Let  t  denote  the  instant 

pa 

of  incipient  plastic  flow  at  the  boundary  x  =  a.  Then  both  t^^  and 
t^^  may  be  determined  frcn  (4.2*5)  and  (4.6.2): 
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since  a  closed  form  solution  for  t  does  not  exist,  ve  nay  deter- 

P* 

mine  t  graphically  by  drawing  a  horizontal  line  of  height  W,  and 
Intercepting  the  temperature  curve  for  the  boundary  x  =  a  as  shown  In 
Fig.  4-9.  Here,  by  (4.6.2), 

w  .  T(a,tj^)  .  (2B  -  -  SM)  Ty(2B  - 

+  2  [2y^  +  B  +  “  2B)]/(2B  - 

For  t>  t  ,  a  new  plastic  Interface,  specified  by  x  «  P,, 

propagates  from  the  boundary  Into  the  Interior  of  the  half -space.  The 

position  Pj  can  be  found  frcni  (4., 6* 2b)  if  we  there  replace  x  by  Py 

amd  t  by  t.  Denoting  the  steady  state  value  of  p,  by  R..  and 
px  p  u. 

setting  T(x,t  )  equal  to  zero  In  (4.6.2),  we  obtain  after  rearrauiglng 
pa 

a  T  E^a  B  +  Et  a/2(l-V, ) 

®11  =  “y^  ■  M  +  E^  a/i-v  ^  (4.6.4) 

We  recall  that  the  steady  state  position  of  the  pleuatically 
deformed  region  is  given  by 

a  T  £  OC 

“r  ■  17 '®  * 

Fran  cui  inspection  of  (4.4.7)  and  (4.6,4),  we  conclude  that  the 
plastic  region  formed  during  unloading  never  extends  to  the  boundary 
R^  of  the  total  plastically  deformed  region. 
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The  total  number  of  regions  that  may  possibly  exist  In  the  half- 

space  as  a  restilt  of  raising  the  naxiotum  boundary  tenqperature  above  . 

the  critlcea  temperature  nay  be  inferred  from  Fig.  4-9 •  The  case 

when  t  is  greater  than  t  _  is  Illustrated  in  Fig.  4-9a.  As  nay 
pa 

be  readily  seen,  in  the  interval  between  t  and  t  only  two  re- 

IU\^ 

glons  exist  in  the  half-spacej  the  elastic  region  of  unloading  and  the 

eleistic  region.  For  t>  t  ,  there  exist  three  regions;  the  region 

pa 

of  plastic  loading  in  tension,  the  elastic  region  of  unloading,  and 
the  elastic  region.  In  Fig.  4-9b  we  have  Illustrated  the  case  when 


t  -  Is  greeter  than  t  .  In  the  interval  between  the  Instants 
mn^^  pa 

tp^  and  t^  ,  four  regions  exist  simultaneously  in  the  hadf -space; 
the  region  of  plastic  loading  in  tension,  the  region  of  elastic  tai- 


loadlng,  the  region  of  plastic  loading  in  compression,  suid  the  elastic 


region.  However,  for  t  >  t  _  only  three  regions  remain,  namely, 
the  region  of  plastic  loading  in  tension,  the  region  of  elastic  toiload- 


Ing,  and  the  elastic  region. 
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Source 

Plane 


Fig.  4-10  Regions  of  Loading  emd  IJhloadlng  In  the 
Plastic -Elastic-Plastic -Elastic  Regime 

4-62  The  Tremslent  Solution  In  the  Region  of  Plastic  Loading  In 
Tension 

We  note  that,  at  the  Instant  t  of  Incipient  yielding  In  the 
region  of  elastic  unloading,  the  sign  of  the  maximum  shearing  stress  Is 
opposite  to  vhat  It  vas  when  yielding  first  occurred  at  that  point,  and 
at  the  Instant  t^.  In  the  present  ceise,  the  stresses  In  the  y  and 
z  directions  will  be  tensile  rather  than  cosq^resslve;  yielding  is  caus¬ 
ed  by  the  continuous  decrease  of  teiig>erature  from  Its  maximum  T8^.ue  T_„. 

iBX 


In  order  to  find  the  stress  emd  strain  In  this  region,  the  equations 
(4.4.20)  and  (4.4.21)  should  be  used  In  the  place  of  (4.4.1^)  and 
(4.4.16).  Recalling 


T  =  M'  T  , 

y  ^ 


(4.4.20)R 


where 


e  =  N*  T  , 
y 


M'  = 


-  2B  (1-b)  -  2E^  g  Up  (4/3) 
(1-s)  +  2  u  (4/3)^^  (1-^j^) 


(4.4.21)R 


(4.4.22)R 


N'  «  [- 


A]  M*  + 


the  stress  and  strain  In  this  region  are  determined  by  Integrating 
(4.4.20)  and  (4.4.21)  between  the  limits  t  and  t^^.  Thus,  by 
(4.6.2),  (4.5.5)  and  (4.5.6),  It  follows  that 

otE,  M'  +  oflB^/l-Z/^  CXE, 

y>,t)  .  M'T(x,t)  .  [(^.  M)  -  a  -  a*)' 

M*  + 

'  ^®m  2B  -  OC&^/l-pj^  *  *^2  y^ 


+ 


(M  -  2B)] 


(4.6.5) 


1+1/j  N'  -  (l+V  /l~P  )  O 

€(x,t)  -  B'T(x,t)  H-  [(»  -  a)  -  a'  -  — 


•(2B  -  -  ac)l  T.  I- 


I  N'  -  (1+Z/^A'^i)  a 
2B  - 


2y, 


N  -  (i+V^/l-V^)  a 
*  B  +  QC&Ji{l-V^ 


(4.6.6) 


The  displacement  U  is  obtained  by  integrating  (4.6.6)  between  the 
limits  X  ewid  I3 


1/2  P, 

U(x,t)  »  a  N'T  (  I  e)  [♦(— ^)  -  »(—■”---)] 


N'  -  (1  +P,/l~'U  )a  OEL 

■*■  ^  2B  -  aE^/1-2^ 


1+2/, 


X  f  N'  -  (1+P  /i-y  )  a 
^  a)]  Tjjj  a  In  -  2B  -  (XE^^/l-Pj^ 


2  y 


B  +  aE^/2(l-^/^) 


N  -  (i+i^^/i-p^)  a 
B  +  0£E^/2(1-Z4) 


(x  “  Pj)  +  g(p3Bt) 


(4.67) 


Here  g(P3B^)  1b  ^  function  of  time,  and  is  determined  by  the  con¬ 
dition  of  continuity  of  displacement  across  the  interface: 
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Aplastic  loeuiing  “  ^elastic  loading  **  *  “  Pj 
In  tension 

(1^.6.8) 


Thus,  by  (4.5.7),  (4.5.10),  (4.6.7)  and  (4.6.8), 


1+Zx, 


a  a  T„(  I  ^ 


2v^ 

1+K 


^  ^  -  1]  +  (N  -  a)  tg  +  0E^/2(1-J'^;  ^^1  “ 


+  a  T  In  -^) 

m  Pg 


a  N  (  I  e) 


1/2 


[♦(- 


2^ 


-♦( 


2V^?^ 


■)) 


(4.6.9) 


for  t  _  >  t  >  t 

pa 


For  the  instants  t  >  t^  ,  the  function  g(Pj,t)  is  the  same  as  that 
given  by  (4.6.9)  except  that  P^  =  Pg  “  Specifically, 


g(p 


1+2/  1/2  p  1+1/ 

3,t)  -5^  aoT^  (  ^e)  [t(-^)  -  1]  +  (n  -  a) 


P, 

2^' 


\  ^  ^  ■  B  +  0®^/2(l-V^; 


(Pj- 


Rj^)] 


for 


t  >  t  > 
pa 


(4.6.10) 
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The  stress,  strain  wd  displacement  for  the  region  of  plastic  load¬ 
ing  in  tension  for  the  cases  t  >t>t^  t_>  t>-t  and 

pa  ^ 

t  >  t  ■>  t  _  are  now  conq?letely  described  by  (h.6.5),  (h.6.6), 
psi 

(4.6.7) ,  (4.6.9)  and  (4.6.10). 

4-65  The  Steady  State  Solution  in  the  Region  of  Plastic  Loading  in 
Tension 

The  residual,  stress  t  auid  e.  are  obtained  by  setting  the 

i\  K 

temperature  T(x,t)  in  (4.6.5)  emd  (4.6.6)  eqml  to  zero,  whereas  the 
permanent  deformation  U_  is  determined  by  letting  the  time  t  in 

A 

(4.6.7)  and  (4.6.10)  tend  to  infinity,  and  by  replacing  by  R^, 

P,  with  R, , .  The  results  are 
5  11 

QE,  M’  +  0®. 

V*>  -  '(« *  r:5i>  -  -  *)1  ■f.  ; 

M'  +  oes  A-K  yjo 

"  2B  -  OE^/l-J/j^  ^o  B  +  aB^/2(l-l/j^; 


(4.6.11) 


e 

r 


(x)  =  [(N 


1+Vj^ 


a) 


N'-  (l+2/j^/L-V^)  a 
2B'  - 


2M)]  T  - 
m  X 


N'  -  {l+V^/l-UA  «  2  y^ 

2B  -  ^o  B  +  0lE^/^l-4y^) 


N  -  (l+vyi-1/3^)  a 

B  +  aE^/^l-Vj_) 


(4.6.12) 
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H*  -  {14-V^/l-l/^)  a 


u^(x)  .  [(N  .  3;^  a)  -  SB  -  aa[yi-z4  . . (2b  - 


N'  -  (1+1//L-2Z)  O  2  y 

■*■  2B  -  B  +  aE^/2(l-i/j^) 


N  -  (l+U^/l-V^)  a 
B  +  aa,/2{i-^> 


>  J-  J- 


y.  W 

+  BT  a£7Fa-%J  ^®i "  ■  r:?5; 


(4.6.13) 


The  stresses,  strains  cuid  displacements  In  the  regions  lying  beyond 
Pj  are  derivable  fixnn  the  relevant  formulas  In  sections  4-3,  4-4  and 

4-5. 


78 


Chapter  An  Elastic  Perfectly  Plastic  Response  of  a  Half-Qpace  to  a 
Uhlformly  Applied  Heat  Pulse 
5-1  Introduction 

The  solution  for  the  response  of  a  half -space  subjected  to  a  uni- 
fojmly  distributed  heat  p^se  was  presented  In  the  preceding  chapter; 
the  medium  was  assumed  there  to  be  elastic,  linearly  strain  hardenlhg, 
and  possessing  a  ten^erature  dependent  yield  stress.  In  this  chapter 
the  same  problem  is  reconsidered  for  the  cases  of  perfectly  plastic  re¬ 
sponse  with  temperature  dependent  yield  stress,  and  also  asstimlng  the 
yield  stress  to  be  constant. 

The  requisite  solutions  may  be  obtained  by  a  direct  calculation, 
or  as  special  cases  of  the  solution  for  strain  hardening  materials. 
Since  the  pattern  of  growth  and  decay  of  various  elastic  euad  plastic 
regions  is  essentlcQly  the  same,  we  shall  not  give  a  detailed  dlscussr 
ion  of  the  solution,  and  limit  the  presentation  to  a  statement  of  the 
solutions  for  stress,  strain  and  displacement  in  the  different  regions. 

The  chapter  is  divided  into  two  parts:  In  Part  1  we  investigate 
an  elastic  perfectly  plastic  medium  with  a  yield  stress  varying  linear¬ 
ly  with  temperature,  whereas  in  Part  2  the  yield  stress  is  taken  to  be 
constant.  The  elastic  solution  is  the  same  for  both  cases,  and  is  siq>- 
plied  by  (4.3.8),  (4  5.9)  and  (4.5.15). 

5-2  Part  1.  Reponse  of  am  Elastic  Perfectly  Plastic  Medium  with  a 

Yields  Stress  Varying  Ljneturly  with  Temperature 
5-21  The  Plastic  Region 

As  before,  we  use  the  yield  criterion  of  Tresca,  so  that  yielding 
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Is  assumed  to  occtu:  vhen  the  maxlnum  shearing  stress  hecomes  eq^jal  to 
the  current  yield  stress  y^  In  shear.  The  -value  of  -the  yield  stress 
y^  Is  new  given  by. 


y^  =  -  BT  (x,t) 

and,  by  (3.3.3)>(3.3.^)  the  yield  functions  assume  -the  fonn, 

^  -  (y^  -  BT) 

where 


. 


Moreo-ver,  since 


T 

y 


2 

have  that 


(5.21) 


(5.22) 


f(l)  .  f(2)  ,  ^(1)  .  p(2) 


(5.23) 


In  -the  present  case  -the  maxlmna  shearing  stress  In  -the  plaatlc 
reglcsi  is  always  equal  to  the  yield  stress;  therefore,  by  (4.3.12), 

Ty(x,t)  =  -  2  [y^  -  Br(x,t)]  (5.2.4) 

The  strains  are  ob-talned  by  (3.3.5)^  Hooke's  Lav,  and  the  condition 
-that  the  toteuL  strains  In  -the  y  eud  z  directions  are  zero.  Sped- 
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fically,  from  (5.2.2)  follow  the  relatlcaw 
s  -  X  =  I  +  I  X^^^  -  I  (X+X) 

Otel  '  X 


(a) 


8  -  8  =  X  X^®^  ^  .  1  x(l)  1  (2)  _  ,  1  ,  /.X 


"y 


^*'1 

E 

«  • 

T  +  a  T 
y 

(c) 

1  -Vi 

E 

•  • 

T  +  a  T 

y 

(d) 

8 

-^y 

(e) 

Moreover,  (5 .2. 4)  inQ>lles  that 


T  «  2  BT 
y 


(f) 


!I^e  solution  of  the  five  equations  (a),  (b),  (c),  (d)  and  (e)  for 
the  five  unkncwna  s^,  and  X  then  yields 

1-2/, 


2  T 


+  2  a  T 


1  -y, 

8  -  X  -  [4B 

X 

*y  “  “I 


+  2  a]  T 


(5.2.5) 

(8) 

(h) 


+  2  a]  T 
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+  a)  T 


1  -K 

=  -  B  =  [2B  —  g"""* 

y  y  ^ 


+  a]  T 


and  the  total  strain  rates  are  e:q)re8sed  "by 


4(1  -  2'U.) 

€  =  s  +£,x  =  [  - = -  B  +  3  Ct  ]T 

X  X 


(5.2.6) 


€  =  0 
y 


(5.2.7) 


The  corresponding  relations  for  regions  in  which  plastic  flow  is 
induced  hy  cooling  can  be  found  without  difficulty,  keeping  in  mind, 
however,  that  dtiring  the  decrease  in  temperature  stresses  in  the  y 
and  z  directions  are  tensile  rather  than  compressive.  In  the  place 
of  (5.2.2)  the  yield  functions  are  now  given  by 


-  (y^  -  ht) 


(y^  -  HT) 


(5.2.8) 


V  -  V 
Z  X 


Since  t  =  t 

y  z 


f(l)  .  f(2)  ,  pd)  .  ,(2) 


(5.2.9) 
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The  stress  and  strain  rates  In  the  region  In  which  plastic  flow  Is 
Induced  by  decrease  In  temperature  may  now  be  expressed  by 


T  =  2  (y  -  HP) 
y 


T  =  -  2  HP 

y 


(5.2.10) 


1  -1/^ 

X  =  [4B  p  '  -  2  a]  T 


(5.2.11) 


4(1  -i/) 

e  =  s  +£.  =  [ - = - 

XXX 


B  +  3  a]  T 


(5.2.12) 


In  order  to  calculate  the  strains  frcm  (5.2.6),  we  Integrate  between 
the  limits  t  and  t^.  Noting  that 


1+V 

e  (x.t,  )  *=  = - rr-  a  T  (x,  t-  ) 

x'  ^  lx'  ^  "^1  ^ 


Is  given  by  the  elastic  solution  (4.5.8),  and  that  by  (4.4.6), 


y 

T(x,tix)  =  T^  = 


(4.4.6)R 


the  total  strain  Is  then  found  to  be 

4(1  -2  2/,) 

r  -L 


»  [- 


E, 


B  +  5  a]  T(x;,t) 


1  +  ly.  4(1  -  2  2/  . ) 

♦ 


B  +  as^i.1  -  V 
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where 

T(x,t)  -  T  e  (5.2.13) 

The  corresponding  displacement  is  obtained  by  integrating  (5. 2. 13) 
between  the  limits  x  emd 

Ml-e^/  )  1/2  p 

u(x,t)  =  [ - +  3a]  a  (  I  e)  [♦(— ^)  -  ♦(— ^] 

h  ®  ^/ct  2  Jet 

1+jy  Mi-2  2/  )  y 

^  B  +  C(Ej_/2(l-  ■  ‘’l^ 

1+Vt  ,  1/2  p, 

+  17^  ®  \  (  2  ‘ 

*^1  “  2  ,/ct 


for  Pj^>  x^a  ,  tjjj>t>tj^ 


or  p^>x>P2 


(5.2.1M 


where  p^  and  Pg  are  determined  from  (4.4.5)  and  (4.5.l)»  re¬ 
spectively. 

5-22  The  gLastlc  Region  of  Iftiloading 

The  stress  and  strain  rates  in  this  region  are  giVim  by 


.  OR 

Ty  (x,t)  -  -  T(x,t)  (4.5.3)R 

i+-y, 

^  T(x,t) 


(M5.4)R 
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Integrating  (4.3*3)  (4.3*4)  between  limits  t  and  t  , 

nx 

where  t>  t  ,  we  find,  by  (4.2.6),  (5*2.4)  and  (5.2.15)»  that 

XfiX 

“  -  Sr  2  vl  - 

(3.2.16) 


1+1/ 


'*<*>*>  '  \  I  * 


Hi-21/  )  y 

+  [  —  B  +  5  a]  -  -  B  +  Cffl^/'2(1-1>^)^ 


(5.2.17) 


Both  (3.2.16)  and  (3*2.17)  are  applicable  for  values  of  x  such  that 

X  >  &,  t  >t>t  and  >  X  >  p,  for  t  >  t  ,  where  t 
2  ^  ^  '  pa  m  2  3  pa'  pa 

and  p^  are  determined  from  (4.6.2a)  and  (4.6.2b)  by  setting  there 
equal  to  zero. 

Integration  of  (3. 2.17)  now  gives 

1+  1/2  Po  P-. 

U(x,t)  =  a  a  T  (  I  e)  [  ♦(~^)  -  *  t(— ^)  -  H 

^  ^1  “  2v^  2.^  2v^ 


1+1/, 


4(1-2  i/)  ... 

+  [ - s-^  B  +  3  a  -  Tp-T-i  a]  [T  a  In  ~ 


®1 


1-2^1  -  '-m  pg 


y„  4(1-2Z/.)  ^  1/2  p^ 

B  +  0®j^/2(l-t/^;  (x-Pl)l  +  [  ®  +  3a]  Tjj  a  (^  e) 


^2'v^. 


(4.2.18) 
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for  P->x>a  ,  t>  t>t 
2  pa  m 

P_  ^  X  ^  P,  ,  t  _■>  t  >  t 
2  3  mR-  pa 


or 


l+'V. 


1/2 


U(x,t)  =  a  a  (  I  e) 


1+7/, 


4(1-21/^)  ^ 

+  [  B  +  3a  -  ®  ^  *  B  +  qE^^2(1-  7/^) 


-  x^)] 


(5.2.19) 


R^:^  X  >  a  , 

t  >  t  >: 
pa  ^ 

^mR. 

R^>  X  >Pj  , 

t>  t  > 
pa 

^mR. 

where  t^,  P^,  R^,  Pg,  in  (4.2.18)  and  (4.2.19)  are  detenninr 
ed  fran  eqwtions  (4.2.4),  (4.4.5)#  (4.4.7)#  (^.5.1)#  (4.5.2)  respect¬ 
ively. 

The  residxial  stress,  strain  and  permanent  defoznation  obtained 
from  (5.2.16),  (5.2.17)  and  (5. 2.19)  a.re  expressed  by 


+  2B)  I  -  2y^  (5.2.20) 

4(1-22/  )  1+1^  y 

ej((x)  =  [  B  +  3a  -  x  "  B  +  aB^/2(l-2/^)5 


(5.2.21) 
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4(1-2  Vj  1+Z^- 

V»)-[— ■■  a] 

Yq  1+  l/2 

+  b;" cffl^72(i-ix^)  (^1-^)1  “ITT];  I®) 

(5.2.22) 

The  equations  (5.2.20),  (5.2.21)  and  (5.2.22)  are  valid  only  In  the 
region  x>  a  for  and  in  the  region  x  ^  R^  for 

Tjj^>  Tjj,  where  is  determined  from  (4.6.3)  by  setting  in  it 
equal  to  zero. 

5-23  Plastic  Reckon  Formed  Durlnp^  Tfaloadlnp 

If  the  maximom  boundary  temperature  is  greater  than  0^, 
given  by  (4.6.3)  for  the  case  when  is  zero,  plastic  flow  occurs 
again  during  unloading.  This  flew  begins  at  the  boundajTr  for  the  in¬ 
stant  t  ,  and,  subsequently,  a  new  elastic -plastic  interface,  sped- 
fled  by  x  =  Py  propagates  from  the  boundary  into  the  interior.  This 
interface  eventxial  ly  reaches  the  position  x  »  R^^^  as  the  teinperature 

T  of  the  medium  decreases  to  zero.  We  let  t  be  the  Instant  at 

px 

which  yielding  occurs  at  the  position  x  during  unloading.  The  values 
■tpx#  tp^,  pj  and  R^  are  given  by  the  following  equations  as 
a  result  of  setting  the  sheeo:  modulus  p  equal  to  zero  in  (4.6.3), 
(4.6.2),  (4.6.2a),  (4.6.2b)  euad  (4.6.4): 

4  y^ 

^  “  (oe^A-^1  +  2b; 


(5.2.23) 


0®. 

(2B  -  T, 


1/2  -X  /V  t. 


a  e 


1 


px  0®. 


PX 


*  fr%*  ; 


4  for  T^> 


V2  -/^V  «, 

V  pa 


=  Jfy, 


for  T  >  T,, 
m  M 


1/2  -  p^l4Kt 

(2B  -  r..  S-S -  e 

1-^1  ^  /w 


*  (3^  t  2B)  I- 

1  5 


4  y  for  T  > 

o  n  w 


a  Tjjj  (c®^A-  +  2B) 

% - ir^; - 

In  the  plastic  region  formed  during  xinloadlng  the  stresses  In 
and  z  directions  are  tensile,  and  e^qpressed  by 

Ty(x^t')  -  T^(x,t)'-  2  ty^  -  HP  (x,t)]  , 
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(5.2.24) 

(5.2.25) 

(5.2.26) 

(5.2.27) 

the  y 

(5.2.28) 


vhereas  the  strain  rate  Is  given  by 
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6^  (x»t) 


[3  a- 


Mi-2  2/J 


B]  T  (x,t) 


(5.?.12)R 


If  the  equation  (5.2.12)  is  integrated  between  the  limits  t  and  t 


px 


where  t>  t  .  then,  by  (5.2.24)  and  (5.2.17),  the  following  result 
px 

is  obtained : 

Mi-2V  )  4(1-2-^/  ) 

e  (x,t)  =  [3  a  - - b3  T(x,t)  +  [ - = - ^  B  +  3  O 

X  El  E^ 


i+t/^  „  y 

^a][T  -  ° 


l-i/^  m  X  ■  B  +  aE^/2(l-l/^) 


]  +  [ 


k{2,r2i/^) 


E, 


B 


1+  P 


» *  irxr  “1  ^  i  A.^/,  [“  -  ■'m  I  *  Sr)l 


for  p,  >  X  ^  a  and  t  >  t 
3  pa 


(5.2.29) 


Integration  of  (5.2.29)  further  gives 


4(1-2  2/^) 


1/2 


U(x,t)  .  [»  -  ■  "  Bl  .  (  I  .n,  -  ♦  (— ^)1 

1  2  2 


4(1-2  2/ J  1+^/.  y 

+  [ — 5 — —  B  +  3a  -  r-rrr  «]  [T  a  m  2 - 8. 


®l 


iZjr  Lijjj  a  xn  -  3  +  cffl^/2(l.  2/^) 


4(1-2  7/)  1+2/  4y 

{x-Pj)^  C  B  -  3a  +  5— a]  taB-oOEj^/i-T.^  ^*"*^3^ 


2B  +  Cffl^/I-T/ 

2B  -  ^  ^  pT  ®  pj^  X  >  a 


(5.2.30) 
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where 


l+y.  ,  1/2  P5  P2  , 

_  a  a  T.  (  5  e)  !♦(— )  -  ♦(— ) 


P  Ml-21/^)  lai'l  »5 

*  ^ ^ ®  *  *“  •  ‘  ^  "pi 

y  Hl-2^J  ,  1/2 

-  B  .  a^/2(i:v-)(-5^i)’-^— ^  ®  *  ‘5 


Pi 

•  t  ♦(— ir)  -  ♦(--Izr)! 

2^/rt  2^ 


for  t 


'mR  >  t  >  t 
1  pa 


(5-2.51) 


Alternately, 


l+V,  ,  1/2  P* 

8(p5.t)  =  iTIT^  a  a  ^  e)  [H—=r)  ‘  ^ 


2v^ 


4(1-2  yj  1+2/,  p, 

+  [ - 5-^  B  +  5  a  -  Ot]  [T„  a  In 


l-y. 


m  El 


B  +  0E^/2(1^Z/'^)  ■  ^5^^ 


for  t  >  t^>  t 


(5.2.32) 


where  t^^^  ,  P2^  Pj  "^pa,  1”  (5.2.29),  (5. 2. 30),  (5.2.31), 

(5.2.32)  are  detennlned  from  (4.2.4),  (4.5*1),  (5.2.26),  (5*2.25), 


respectively. 
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The  residual  stress,  strain  and  permanent  deformation  for  this 
region  are  found  from  (5*2,28),  (5*2.29),  (5*2*30)  and  (5*2.32)*  Thus 
we  obtain  the  following  results 


=  2  y 
R  •'o 


(5*2.33) 


1+2/, 


4(1-2  ly.)  a  Jq 

ep(x)  =  [  g-  B  +  3a  -  X  "  B  +  cffl^/2(l-  ^ 


+  [ 


4(1-2  2y)  1+-2/. 

•A- 


B  -  3a  + 


l-p^  2B  -  E^a/l-y^ 


OE^ 

T  -  (2B  +  ri;7")3 


(5*2.34) 


4(1-21/.)  1+2/. 

Vx)  =  [-^B+3a-^  a]  [\aln- 


y  4(1-2  U  ) 

B  +  aE^/2(l-Z/^)  B  -  3a 

1+2/  2B  +  aE/l-K  R„ 

*-2B  -  aE^/l-2/3^  ®  ^  T" 


2B  -  OEj^/l-  ^\l  " 


(5*2.35) 


Equations  (5*2.25),  (5*2*54)  and  (5*2.35)  are  valid  only  for  the  region 
R^j^>  X  >  a  where  R^j^  Is  determined  hy  (5*2.28). 
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5-3  Part  2.  Response  of  an  Elastic  Perfectly  Plastic  Medium  with  a 

Constant  Yield  Stress 

If  the  infinite  half -space  is  subjected  to  a  uniformly  applied  heat 
pulse  at  its  boundary,  and  if  the  medium  is  assvnned  to  be  elastic,  per¬ 
fectly  plastic  with  a  constant  yield  stress,  the  solution  is  readily  ob¬ 
tained  as  a  special  case  of  the  results  described  in  Part  1  of  this 
chapter.  Specifically,  the  equations  describing  the  stress,  strain  and 
displacement  in  various  regions  for  the  medivun  considered  here  may  be 
deduced  by  setting  the  quantity  B  equal  to  zero  in  the  corresponding 
equations  derived  in  Part  1  of  this  chapter.  The  resrjlts  of  this  sim¬ 
plification  will  be  presented  in  the  subsequent  sections,  the  elastic 
solution  being  again  given  by  (k.^.8),  (4.3*9)  and  (4.3*15) • 

5-31  The  Plastic  Region 

In  this  region,  by  (5*2*10),  (5*2.13),  (5*2. l4),  the  stress, 
strain  and  displacement  are  expressed  as  follows: 


Ty(x,t)  =  -  2  y^ 


1+2/  2(1-2/) 

^(x,t)  =  30!  T(x,t)  -  (3  -  jrjT^)  - 


E,  ^o 


(5*3*1) 

(5*3*2) 


C(x,t)  =  a  a  T  (  i  e)  [3  [t( - - — )  -  If  ( - ^ — )] 

“  ^  ^  2  ^t  2  ^ 


1+ 1/,  P, 


1*2/.  2(1- t'.) 

*  (“1  - 

■^"1 


for  P3^^x>a  ,  t^^t^t^ 


P,  >  X  ^  p„  ,  t_„  >  t  t 


mR, 


m 


(5*3*3) 
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where  Pg  is  determined  from  (IJ..5.I)  and  from  (4.4.5)  t>y  setting  B 
there  equal  to  zero. 

5-32  The  Elastic  Region  of  Unloading 

From  (5.2.16),  (5.2.17),  (5.2.18)  and  (5. 2.19)  follows  that  the 
stress,  strain  and  displacement  in  the  elastic  region  of  unloading  may¬ 


be  written  as 


ctJS,  _ 

T  (x,t)  =  -  Y^pr-  [T(x,t)  -  2y^  (5.3.4) 


1+U  1+2/,  ^  2ii-y) 

(5  a  -  5^  Qt)  [T^  -  -  yj 


(5.3.5) 


Both  (5.3.4)  emd  (5.3*5)  are  applicable  under  the  conditions 


p-  >  X  >  a  ,  t  „  >  t  >  t 


Pg  ^  X  >  Pj  ,  t  >  tp^ 


u(x,t)  =  Ypr  ®  O'  ("  71:)  - 


+  (3  Of  -  1-7^0')  I  ®  \  [>^(— ~)  -  (— ~)] 

i  ^  2  gy;j^  gyij^ 


^  2(144 


a  yo  (Pl  -  x)  "  a  in  - 


Pp  >  *  >  >  % 


(5.3.6) 
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or 


1+K  1/2 

U(x,t)  -  a  a  T  (  I  e)  [♦(— ^)  -  1] 

1+1/  2(1-1/) 

+  (3  a  -^a)[  (R^  -x)  +  a  In  —  ] 


for 


,  t  >t>t^ 
1  ^  pa  mR 


R,>  x>P,  ,  t>t  >  t  „ 

1  3  pa  mR^ 


(5.5.7) 


where  P^,  P^,  ,  R^,  and  P^  in  (5.3.^)  to 

(5.3.7)  are  obtained  from  (i4^.4.5)>  (^.5.l)>  (^.2.4),  (4.5*2), 

(4.4.7) ,  (5. 2'. 25)  and  (5.2.26),  respectively,  with  the  quantity  B  in 
these  eqmtions  being  set  equal  to  zero. 

From  (5.5.4),  (5.3.5)  and  (5.5»7)  've  now  obtain  the  following 
results  for  the  residxial  stress,  strain  and  displacanent : 


2(1-Z^, ) 

SW  -  (»  Wm?- 

E^a 

^o^ 

l.K 

2{l-1{) 

U„(x) 

'  1''.“ 

X 

E^a 

(R^-x 

(5.3.8) 

(5.3.9) 


l+i^l  „  1/2 

WT"”.'!*) 


(5.5.10) 


Equations  (5*3.8),  (5*5.9)  and  (5.5*10)  are  valid  in  the  region 
R^:^  X  >  a  for  ,  or  in  the  region  R^^  x  ^  R^^  for 

Tj^>  where  by  (5*2.27),  (5.2.25),  and  Tj^  are  expressed  in  the 

present  case  by 


9^ 


a  T 


'm 


ill 

o«i 


(5.3.11) 

(5.3.12) 


5-33  The  Plastic  Region  Formed  During  Iftiloadlng 

If  the  majclmuni  houndary  temperature  is  greater  than  !I^  given 
by  (5.3.12),  a  new  plastic  region  of  unloading  emerges  from  the  boundary 
at  the  instant  t  .  This  region  eventually  extends  to  the  position 
X  =  which  is  given  by  (5.3.11).  The  relevant  solutions  for  this 
region  are  found  from  (5.2.28),  (5.2.29),  (5.2.30),  (5.2.21),  (5.2.52), 
and  are  as  follows: 

Ty(x3^t)  -  2  y^  (5.3.13) 


1+P 

e3j(x»t)  »  3  a  T  (x,t)  +  (5  a  -  y^ 

U(x,t)  =  3  a  T  a  (  I  e)  (♦(— ^)  -  t(  ^  __)) 

“  2^  2  ^ 

1+p  2(1-;/) 

+  (5a  -  j^jroc)  -  y^  (x  -  p^)  +  g  (P5,t) 


(5.3.U) 


for 

where 


P,  >  X  ^  a 
3 


1+K 


1/2 


(5.3.15) 


s<V>  =  rrjf 

1+^1  f  1/2  Pp  p 

2(1-;/ )  p  -1 

+  — sr77“  y  (Pi  ~  Pi)  +  T  a  Ip  «  '  for  ^  r,^  t  >  t 

E^a  'o'!  3  m  Pg  j  mR, 


pa 

(5.3.16) 
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^  (  I  -  1)  .  (3  a  -  ^a) 

P  2{1-V) 

[T^a  ln-2+  (R^  .  p^)] 


for 


t  >  t  >  t  ^ 
pa  OR, 


(5.5.17) 


It  should  be  noted  that  p,  and  t  in  (5. 5.16)  and  (5.5.1?) 

p  pa 

can  be  determined  from  (5 •2,26)  and  (5.2.25)  with  the  qtuuatity  B  there 
being  set  equal  to  zero.  Finally,  the  residual  stress,  strain  and 
permanent  deformation  obtained  from  (5.5.15),  (5.5.1^),  (5*5.15)  and 

(5.5.17): 


-  2  y^ 


(5.5.18) 


•L+V  2(1-2/) 

^W-(3«-r:^«)-^yo  (5.5.19) 

1.+U  R  , 

l^(x)  =  (5a-jr^a)  T^aing^ 

2(1-7/)  l+Z^  1/2 

^  E^a  y©  (*  +  "  2Rii)  "  \  ( 2 

(5.5.20) 

The  equations  (5. 5*18),  (5*5.19)  and  (5.5.?0)  are  applicable  only  in  the 
region  x  ^  a  where  R^ 


is  given  by  (5. 5*11). 
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(Jhapter  6.  Elastic  Plastic  Besponse  of  a  Laterally  Constrained  Plate 
to  a  Italformly  Applied  Heat  Pulse. 

6-1  Introduction 

Transient  thermal  stress  analysis  of  free  plates  has  "been  treated  hy 
Weiner,  Landau,  Zwlcky  [6.1]  [6.2]  [6.5],  Yushel  [6.4],  sind  most  recent¬ 
ly  by  Mendelson  and  Spero  [6.5],  In  this  chapter,  we  shall  consider  the 
response  of  an  Infinite  plate  constrained  In  the  lateral  direction,  and 
subjected  to  a  xmlformly  applied  heat  pulse  at  one  boundary,  the  other 
boundary  being  held  fixed  and  either  maintained  at  zero  temperature  or 
Insulated. 

The  ptirpose  of  studying  a  plate  of  finite  thickness  is  to  remove  a 
basic  shortcoming  of  the  solution  found  for  the  half -space,  namely,  the 
absence  of  a  characteristic  dimension  of  length.  SpeclflctQly,  we  shall 
relate  the  dimension  of  the  plasticedly  deformed  region  to  the  plate 
thickness  and  thus  establish  a  basis  for  obtaining  approximate  solutions 
of  plate  problems  from  the  corresponding  solution  for  the  half -space. 

If  the  mechanical  and  thermal  boundary  conditions  for  a  plate 
problem  are  similar  to  those  applied  for  the  h£J.f -space,  then  the 
analysis  will  also  be,  in  essence,  identical  with  that  carried  out  for 
the  half -space.  In  fact,  the  specific  solution  derived  in  Chapters  4 
and  5  also  apply  directly  for  plates,  provided  that  the  tenQ)erat\xre 
distribution  for  the  half -space  is  replaced  by  the  appropriate  tempera- 
tTire  distribution  in  the  plate.  Therefore,  let  us  first  investigate  the 
tenperature  solutions  for  plate  problem. 

6-2  The  Temperature  Problems 

The  tenperature  distribution  in  an  infinite  plate,  produced  by  a 
heat  pulse  applied  uniformly  at  one  boxmdary,  the  other  boindary  being 


maintained  at  zero  temperature^  may  be  obtained  by  sviperposlng  the 
teniperature  fields  induced  by  a  suitably  located  instantaneous  source 
and  sink. 
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j  Plane  Sink 
I  or  Source 


Pig.  6-1.  Location  of  Source  and  Sink  for  a  Plate  Subjected  to  a  Uni¬ 
formly  Applied  Heat  Pulse  at  One  Boundary,  With  the  Other 
Boundary  Either  Maintained  at  Zero  Temperature  or  Insulated. 

Let  the  plate  be  initieOly  at  a  uniform  temperature  T^  .  We 
characterize  the  origin  of  coordinates  as  shewn  in  Pig.  6-1,  and  assrone 
that  the  boundary  x  =  o  is  fixed  while  the  other  boundary  at  x  »  L 
remains  free.  An  Instantaneous  plane  source  and  sink  of  equal  strength 
are  placed  symmetrically  to  the  x  =  o  plane;  due  to  the  symetry  of 
these  locations,  the  teniperature  change  at  x  s  o  is  zero.  At  the  free 
boundary,  or  at  emy  other  plane  between  the  two  boundaries,  the  tenpera- 
ture  will  Increase  from  zero  to  a  maximum  emd  again  decrease  to  zero. 
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The  8TQ>erpoBitlon  of  temperature  solutions  due  to  a  source  and  a  sink  to 
describe  the  Increase  of  tenqperattire  over  the  reference  temperature 
Is  thTis  expiressed  by 


T(x,t)  «  C. 


1 


[e 


(7lrf-L-x)‘ 
■  Wt 


-  e 


Wt  ] 


(6.2.1) 


Here  Is  a  constant,  L  denotes  the  thickness  of  the  plate,  and  yL 
represents  the  distance  of  the  source  from  the  free  boundary  (Fig.  6-1), 
the  quantity  7  being  an  arbitrary  constant.  By  varying  7  >  the  shape 

of  the  heat  pulse,  or,  equlvsQ-ently,  the  time  rate  of  Increase  of  tentper- 
atiire  may  be  adjtosted.  By  Inspection  of  (6.2.1)  we  find  that  T  «  o  for 
X  ■«  o.  If  two  sources  of  equal  strength  are  placed  sjnmmetrlcally  to  the 
X  a  o  plane,  then  no  heat  will  be  transferred  across  this  plane,  and  the 
ten5>erat\ire  distribution  will  correspond  to  the  case  when  the  plate  Is 
subjected  to  a  uniformly  applied  heat  pulse  at  the  free  boundary  •(rtiile 
the  fixed  boundary  is  Insulated  (Pig.  6-I).  We  find 

(yl/fL-x)^  _  (7]>]>x)^ 

T(x,t)  -  Cg  [e"  +  e  "  ^ttct  ]  (6.2.2) 

t 


where  is  a  constant. 

The  maximum  temperature  T  attainable  at  an  arbitrary  plane  x 

mx 

in  the  plate,  and  the  corresponding  instant  t^^^  at  which  the  tempera¬ 
ture  assumes  this  value  T  are  very  significant  in  the  analysis  of  the 

mx 

mechanlceil  response.  In  order  to  evaluate  T^jj^  ,  it  is  first  necessary 
to  determine  t  frcm  the  condition 


dr\ 


o 


(6.2.3) 
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The  solutlona  appropriate  to  (6.2.1)  and  (6.2.2),  respectively,  are 
then  given  by 


[e 


mx  Sk  mx  2k 


mx 


_  (yl^h-x)^  _  (rLfL-x)^ 

e  ”  Wt  -  e  "  iwct 

mx  mx 


for 


L  >  X  >  0 


(6.2.4) 


and 


_  [e 


mx 


mx 


2k 


_  (yi^b-x) 

"TJJft  +  e 

mx  mx 


(yl>]>x)‘ 

Iwct 


for 


L  ^  X  ^  o 


(6.2.5) 


The  constants  and  In  (6.2.1)  and  (6.2.2)  may  be  eliminated 
by  introducing  the  maximum  values  T^^  ,  ordinarily  assumed  to  be  known, 

for  the  plane  x  *»  L  .  Denoting  by  t^  the  Instant  for  which 
T  =  ,  we  obtain  from  (6.2.1),  (6.2.4),  (6.2.2),  (6.2.5)  the  results 


[e 


.m 


mL 


i2l£ 

2K 


,  ,2  , 

fe23iLii/2 


[e 


]V2 


(6.2.6) 


ICX) 


where  in  (6.2.6)  and  (6.2.7)  is  found  from  (6.2.4)  wd  (6.2.5)# 

respectively. 

Having  derived  expressions  for  and  Cg  ,  we  are  now  in 

a  position  to  calculate  the  msocinnm  temperature  attainable  at  any 

point  in  the  plate.  The  variations  of  temperattire  T  with  time  t  for 

various  values  of  x  for  both  solutions  (6.2.1)  and  (6.2.2)  have  been 

plotted  in  terms  of  dimensionless  variables  t/p^^  #  x/l  ,  and 

2 

t/ii  /2x  ,  and  are  shown  in  Figs.  6'»2  and  6-5#  respectively. 

6-3  Response  of  a  Laterally  Constrained  Plate  to  a  Uniformly  Allied 


FIG  6-2  TEMPERATURE  CHANGE  WITH  TIME  AT  DIFFERENT  POSITIONS 
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FIG  6-3  TEMPERATURE  CHANGE  WITH  TIME  AT  DIFFERENT  POSITIONS 
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The  Blmllarlty  between  the  problem  for  &  plate  auxd  the  half -space 
becomes  evident  from  Inspection  of  equations  (4.3«5)  to  (4.5.7).  In 
peurtlcular,  comparison  of  the  temperat\xre  distribution  r^resented  by 
Pig.  4-2  and  Pig.  6-2,  and  the  fact  that  the  temperature  and  displacement 
In  the  half -space  were  assumed  to  vanish  for  large  x  ,  bring  out 
clearly  the  analogy  between  the  two  problems.  Therefore  the  response  of 
the  plate  to  the  mlfoimly  applied  heat  pulse  represented  by  (6.2.1)  is, 
In  essence,  the  same  as  that  of  the  half -space.  Especially,  If  the 
proper  tenqperature  functions  are  used,  the  equations  describing  the 
transient  and  residual  stresses  and  strains  In  the  half -space  may  be  used 
to  describe  the  corresponding  qtiantltles  In  the  plate. 

During  the  Initial  stage  of  the  heat  pulse,  the  entire  plate  Is 
elastic.  The  non  vanishing  stress  and  strain  components  are  given  by 

(4.3«9)  «uid  (4.5.8): 


T 


T  =  T  = 

yy  zz 


a  £L 


e 


e 

XX 


a  T 


(4.5.9)  Re 
(4.3.8)  Re 


The  maximum  shearing  stress  Is  represented  by 

T  °®1 

^  °  2'(l-j^)  ^  (4.3.12)  Re 

irtiere  the  temperature  T  for  (4.5.9)  Re,  (4.3.8)  Re  and  (4.3.12)  Is 
given  by  (6.2.1). 

The  displacement  Is  obtained  by  Integrating  (4.3.8)  Re  In 
c<HiJ\uictlon  with  (6.2.1).  In  view  of  the  condition  (6.3.4)  we  obtain  the 
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result 


U. 


(x^t) 


a  (»c)V2  [2*(J2^  .  ♦(24^  .  t  (jaa^ )] 

2^  2y/Kt  2y/!a 

(6.3.5) 


Source 

Plane 


Source  of 
Sink  Plane 


Fig.  6-4  BoxuidarJes  of  the  Total  Plastically  Deformed  and 
Steady  State  Plastic  Regions  in  the  Plate. 

It  is  seen  from  Figs.  6-1  and  6-4  that  the  origin  of  coordinates  is 
chosen  at  the  fixed  boundaiy  of  the  plate  instead  of  at  the  source 
plane  as  was  done  in  the  case  of  the  half -space  (Fig.  4-10).  Since  the 
heat  pulse  is  applied  at  the  plane  x  =  L,  the  dimension  of  the  toted 
plastically  deformed  region  and  that  of  the  steady  state  plastic  region 
are  given  by  (L  -  R^^)  and  (L  -  respectively.  Here  euid 

still  denote  the  distances  between  the  origin  end  the  boundary  of 
the  corresponding  regions.  In  order  to  estimate  the  damage  produced  by 
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the  heat  pulse,  the  location  of  the  boundary  of  these  two  regions  will 
now  be  determined.  By  the  criterion  of  Tresca,  yielding  occurs  where - 
ever  the  maximum  shearing  stress  becomes  equal  to  the  yield  stress  in 
sin5)le  shear.  From  examination  of  (4.3.12),  (6.2.1)  and  Fig.  6-2,  it 
is  evident  that  the  Initial  yielding  starts  at  the  boundary  x  =  L 
(Fig.  6-4),  and  that  the  elastic  plastic  interface  will  progress  into 
the  plate  from  this  boundary  as  the  temperature  continues  to  Increase. 
Following  the  same  reasoning  as  that  presented  in  section  4-4,  the 
boundary  of  the  total  plastically  deformed  region  is  thus  determined 
by  solving  the  following  two  simxxltaneous  algebraic  eqxiations : 


7LtL-R. 


_  (7L*-LfR^)‘ 


mR 


[e 


1  (7L»-L-Rj^)‘ 

ac 


4irt 


mR. 


-  e 
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2k 


mR 


^  _  (7L+L-R^)^ 


_  (71+I/<-Rj,)‘' 


'mR, 


-  e 


TJct 


(6.5.6) 


mR, 


emd 


B  +  E^a/2(l-i/^)  ■  ^mR. 


(6.3.7) 


where 


(7]>L-R^)  -  (7LfI/fR^)' 


Iwct 


T  =  C  1  ,  “^1 

irR^  '^l  7  1/2  le 

^  ^mR, 


4xt 


niR, 


Here,  it  may  be  noted  that  (6.3.7)  Is  essentieJ.ly  identical  to  (4.4.6). 

Simple  and  explicit  expression  for  R^  in  terms  of  other  known 
quantities  cemnot  be  found  from  (6.3.6)  and  (6.3.7)i  therefore,  the 
following  graphicail  method  is  suggested  (Fig.  6-5). 


io6 


Fig.  6-5 


Graphical  Determination  of  R^, 


andR^l 


By  the  use  of  (6.2.1)  and  (6. 2 A),  the  corresponding  values  of  t^ 

and  T  for  a  given  value  of  x  may  he  calculated,  and  plotted  as 
mx 

shown  in  Fig.  b-5.  If  a  horizontal  of  height  equal  to  T  _  and  given 

by  (6.3.7)  is  drawn,  then  the  coordinate  of  its  intersection  with  the 

T  curve  gives  both  t  „  euid  R^ .  In  the  case  when  the  maximum 
mx  mR^  1 

boundary  ten^jerature  is  greater  than  for  which  plastic  flow 

occurs  during  unloading,  a  new  plastic  region  would  be  fomed  in  the 
region  of  elastic  -unloading  and  the  steady  state  position  R^^^  of  this 
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plastic  elastic  interface  may  be  deteimined  by  locating  in  the  plasti¬ 
cally  deformed  region  the  plane  at  which  the  maximum  tejjg>erature  attain¬ 
able  is  equal  to  T,- .  We  obtfiin  the  solutions  for  the  plate  in  the 
plastic  region  formed  during  unloading  by  replacing  in  {k.6.^),  (4.6.6) 
the  tenq)erature  solution  appropraite  to  the  half -space  problem  by 
(6.2.1).  In  particular,  replacing  the  mEocimm  tenqierature  a/x  by 
in  (4.6.2),  we  find  that 

^  W,  ^o  (M-2B)] 

M  +  — = - B  B  +  — = - 

2(l-p^)  2(1-2/^)  (6.5.8) 


By  drawing  in  Fig.  6-5  a  horizontal  of  height  equeLL  to  given  by 

(6.5.8),  the  value  of  can  be  read  off  directly  for  any  value  of 

7  from  the  point  of  intersection  of  this  horizonted  with  the  T 

mx 

curve.  This  process  is  equivalent  to  the  simultaneous  solution  of  the 
following  two  equations  for  and  t^  : 


(7LtL-R^^)' 


(7lH-LfR^^)' 


Tpft 


'I  t 


[e 


mR 


T3ct 


11 


mR. 


-  e 


11 


’mR. 


11 


oos. 


OE., 


(M-2B)] 


M  + 


-  B 


2(1-^/^) 


B  + 


2(1- 


(6.5.9) 
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”®11  (21/4-L-R^^)^  ”®11  (71/Hh-R^^)^ 

^  _  2k  ~  ~  e _  2/C  ^ 

”^11  _  (rlri-L-R^^)^  _  (ylH-LfR^^)^ 

Wt  „  Jwct  _ 

®  “®ii  "  ^  ”®ll  (6.5.10) 


The  quantity  in  (*t-.3.9)  is  given  by  (6.2.6).  With  R^  and  R^^ 
determined,  the  position  of  the  steady  state  elastic  plastic  interface 
and  the  dimension  of  the  total,  plastically  deformed  region  are  known. 
The  remaining  steps  in  the  plate  problem  are  entirely  analogous  to 
those  of  the  half -space  problem,  and,  therefore,  will  not  be  elaborat¬ 
ed  here. 

6-4  The  Response  of  a  Laterally  Constrained  Plate  to  a  Heat  Pulse 
Uniformly  Applied  at  its  Free  Boundary  with  the  Fixed  Boundary 
Insiilated 

The  elastic  solution  for  the  present  case  is  the  same  as  that  pre¬ 
sented  in  the  last  section,  except  that  the  temperature  function  T  in 
(4.3.9)>  (4.3.8)  and  (4.5.12)  is  now  given  by  (6.2.2).  The  counteiTpart 
of  (6.3.5)  then  reads 

U(x,t)  =  rt  c  («/t)^'^^  t(23i!i^)] 

1  2./cb  2yKt 

for  X  ^  0  (6.4.1) 

If  the  temiperature  T  „  given  by  (6.3.7)  is  greater  than  T  ,  which 

mo 

is  defined  as  the  maximum  temperature  attainable  at  the  fixed  and  insu- 
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lated  boundary,  the  position  of  the  toted,  plastically  deforaed  region 
may  be  obtained  from  the  solution  of  the  simultaneous  equations 


(7LfL-R^)‘ 


Wt 


(rLfLfR^)' 


mR, 


iwct 


mR., 


[e 


(6.4.2) 


and 


^(7L+L-R^)‘ 
Wt 


_(7lH-LfR^)' 


”®1  (7lrfL-R3^)^ 

S'  "■'a, 


Wt 


mR. 


1  (7L+IH-R, ) 

- 


1  _  (7]>L-R^)‘^  _  (7l/4-L+Rj^)‘ 

^  ®  '  -TScT - 


mR, 


mR, 


(6.4.3) 


where  in  (6.4.2)  is  » iven  by  (6.2.7).  The  graphical  method  sxig- 

gested  in  the  last  section  for  locating  R^  is  still  applicable  here. 

If  on  the  other  heuid,  T  _  is  less  theui  T  ,  the  entire  plate  would 

'  mRj^  mo'  ^ 

e^erience  varying  amounts  of  plastic  loBuilng  and  would  became  elastic 
agedn  once  the  temperature  begins  to  decrease. 

The  maximum  boxindary  temperature  !I^  for  the  present  case  is 
also  given  by  (6.3.8). 

The  steady  state  position  R^^  of  the  elastic  plastic  interface 
in  the  plate  for  meucimum  boxmdary  tenq>erature  Sweater  than 

suid  for  a  given  value  of  7  is  obtained  from  the  solution  of  the 


no 


simultaneous  equations 


Iwct  ^  ”  hct 


loR 


mR 


T72 


[e 


11 


mR 


+  e 


'll 
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OE, 


[2  + 


(M-2B)] 


M  + 


2(l-^l) 


B  + 


2(1-  2/^) 


(6.4.4) 


and 


.(/LfL-R^^) 


_  (rUI/fR^) 


iMCt 


mR, 


^11  (rlrf-L-B^J 

t»  --IT-—- 


4ct 


mR 


11  (ylH-LfR^)' 


2k 


"®n  _  (ri>L-R^)^  _  (rLfUR^) 


iwct 


+  e 


mR 


Iwct 
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mR 


n 


(b.4.5) 


where  In  (6.4.4)  Is  given  by  (6.2.7).  The  value  of  R^  can  easo 
be  determined  graphically  a.8  was  Indicated  In  the  last  section. 

The  transient  and  residual  solutions  for  the  plate  considered  here 
can  be  derived  In  the  same  manner  as  was  done  In  Chapters  4  and  No 
additional  chanenge  Is  posed  by  these  problems,  hence,  further  analysis 
Is  not  pursued. 


Chapter  7  Ihmerlcal  Results  and  Dlscusslnn 
7-1  Introduction  and  General  Discussion 


m 


A  shortcoming  of  the  half -space  problem  is  that  it  lacks  a  charac¬ 
teristic  dimension  of  length;  In  order  to  remove  this  shortcoming  a 
plate  of  finite  thickness  was  treated  in  Chapter  6.  We  noV  consider  the 
question  of  \rtiether  it  is  possible  to  predict  the  transient  and  residual 
stresses  and  deformations  in  a  plate  from  the  solution  of  the  half -space 
problem.  This  question  is  considered  in  some  detail  in  the  following 
pages. 

The  analogy  of  problems  for  plates  euid  the  half -space  was  previous¬ 
ly  referred  to  in  Chapter  6,  and  it  was  noted  there  that  the  equations 
for  stresses  and  strains  in  the  half -space  problem  may  be  used  for  the 
plate  problems  if  the  tetnperature  function  is  appropriately  modified. 
Therefore,  it  is  of  interest  to  determine  the  extent  to  which  special 
results  of  the  half -space  problem  may  be  directly  applied  to  plate 
problems . 

In  the  case  that  the  effect  of  the  dependence  of  yield  stress  on 
tenqperature  dominates  over  the  strain  hardening  effect  of  the  material, 
the  maximum  transient  stress  attainable  in  the  half -space  is  found  from 
(4.5*8),  (4.4.4),  (4.4.6)  to  be  dependent  only  on  material  properties, 
and  is  equal  to  twice  the  initial  yield  stress  in  shear.  The  maximum 
trauislent  stress  for  the  plate  is  therefore  given  by 


where 


(7*1.1) 
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A  close  inspection  of  (4.6.11)  describing  residual 

stress,  and  of  (4.4.27),  (4. 5-15),  (4.6.11)  describing  transient  and 

residual  strains,  reveals  that  these  quantities  depend  in  addition  to 

material  properties,  solely  on  the  maximum  temperature  attainable  at  a 

particular  position,  and  assume  their  maximum  values  at  the  boimdary 

X  =  a.  The  highest  temperature  T^  attainable  in  the  haLlf -space,  emd 

the  highest  tenqjerature  attainable  in  the  plate  are  ordinarily 

assumed  to  be  known.  If  T  is  set  equal  to  T  the  maximum  trems- 

m  tni. 

lent  strain,  maxlmtim  residual  stress  and  strain  of  the  half -space  be¬ 
come  equal  to  the  corresponding  quantities  of  the  plate.  Therefore,  in 
view  of  the  foregoing  remarks,  and  by  (4.4.27),  the  maximum  transient 
strain  in  the  plate  may  be  expressed  in  the  form 

1+y,  y 

®max  =  N  Tjjj  -  [(N  -  B  +  a  Ej^/2(1- //^) ^  (7-1-2) 


whereas,  by  (4.5.12),  (4.5-15),  the  maximum  residual  stress  and  strain 
are  obtained  from 


max 


y 

=  '^mL"  ^  B  +  0tE^/(2(l-i/^) 


=  (N 


1+1/, 


max 


(7.1.3) 


T  T 
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Here  is  the  maxlnjum  free  boundary  temperature  for  which  plastic  flow 
in  the  plate  during  unloading  will  impend.  For  the  case  when  1^, 

it  follows  from  (4.. 6. 11),  (4.6.12)  that 

a  E.  M*  +  a  a  E, 

'«max  “  ^  ’^”1' 


M'  +  a  K./1-U  y 

■  2B  -  a  ^^o  ■*■  B  +  a  E^/2(l-/^) 


^Wx 


1+P^  N‘  -  1+P/l-U  a  E- 

■  2B  -  a  E^A-^i  ’  1^  ” 


(7.1.4) 


*jnL 


f"'  -  (i.^iA-H)  °  [i.  V  . 

2B  -  a  Ej^/l-2£ 


2  y. 


^o  B  +  a  Ej^/2(l-  ^) 


N  -  {l*P^/x-P^  a 
B  +  a  E^/2(l-^)  ^o 


It  is  thus  seen  that  the  maximum  transient  and  residual  stresses 
and  strain  induced  in  a  plate  by  a  heat  pulse  associated  with  the  maxi¬ 
mum  boundary  ten^perature  may  be  predicted  by  conQ)uting  the  corres¬ 
ponding  quantities  for  the  half -space  subjected  to  a  pulse  of  maximum 
value  T^.  We  also  note  that  the  shapes  of  the  two  pulses  need  not  be 
the  same,  but  that  their  meucimum  amplitude  must  be  equal  in  order  for 
the  similarity  to  exist. 

Equations  (7.1.1)  to  (7.1.4)  together  with  (6.3.6),  (6.3.7)# 
(6.3.9)#  (6.3.10),  (6.4.2),  (6.4.3)#  (6.4.4)#  (6.4.5)  for  the  determi¬ 
nation  of  the  positions  and  of  the  boundary  of  the  total 
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plastically  defomed  and  steady  state  plastic  regions  may  te  consider¬ 
ed  to  r^resent  the  significant  parts  of  the  transient  and  steady  state 
solutions  for  the  plate  problems. 

It  is  also  possible  to  obtain  the  residual  stress  and  strain  at 
any  other  place  in  the  plastically  deformed  region  of  a  plate  by  re¬ 
placing  T  -  in  (7.1.3)  or  (7.1.4)  by  the  corresponding  maximum  tem- 
peratxire  attainable  at  that  position.  Since  the  dispiacemenT, 

boundary  conditions  for  the  two  problems  are  not  exactly  identical, 
the  solutions  for  displacement  will  differ,  in  essence,  by  a  rigid 
translation. 

We  adopt  the  qxaantities 

TAni»  PiA>  PjA^  RiA>  ®iiA;  ^A  «  t/a^/2/( 

to  represent  dimensionless  temperature  T,  position  x,  interface 
positions  p^,  Py  R^j^,  displacement  U  and  time  instant 

t,  respectively,  for  the  half -spe.ee  problem.  Here  a  is  the  distance 
between  the  source  plane  and  the  boundary  of  the  half -space;  for  con¬ 
venience,  the  dimensionless  time  is  often  denoted  by  n. 

Similarly,  the  quantities 

^A'y  ^llA 

are  adopted  to  represent  the  dimensionless  temperature  T,  position  x 
and  steady  state  interface  positions  R^  and  R^  for  the  plate  prob¬ 
lems.  The  dimensionless  stress  and  strain  for  both  cases  «a*e  expressed 
by  t/y^  and  c/fafT,  respectively. 
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A  detailed  calcvilatlon  of  stresses  emd  deformations  In  a  half- 
space,  and  for  specific  data,  vlU  enable  us  to  gain  an  Insight  Into 
the  general  nature  of  the  thermcmechanical  response.  In  view  of  the  es¬ 
tablished  slmlLarlty  with  plate  problems,  we  shall  also  be  In  a  position 
to  draw  concliision  regarding  the  behavior  of  plates. 

An  almlnum  alloy  (25ST)  having  the  following  material  properties 
was  selected  for  numerical  calculation: 


=  10  X  10^  psi  , 

K  =■  0.133  in^/sec.  , 

P  =  0.35 
P 

=  3.9215  X  10^  psi, 

y  =  -  BT 

where  y^  =  23,000  psi  , 


a  =  14.5  X  10'^  ln/in/^°  , 

2^1=  0.275  , 

Up  =  0.2891  X  10^  psi  , 

B  =  0.293  X  10^  psi/i*®. 


(7.1.5) 


The  data  (7.1.5)  repi'esent  the  thenmd  and  mechsuiical  properties  of  an 
elastic,  linearly  strain  hardening  medium  having  a  yield  stress  in  shear 
which  varies  linearly  with  temperature. 

In  order  to  compare  the  solutions  presented  in  chapters  4  and  5,  a 
perfectly  plastic  medium,  with  thennal  auid  mechanical  properties  identi¬ 
cal  to  tnose  given  in  (7-1.5)  was  also  considered,  except  that  there  the 
plastic  shear  .modulus  and  the  coefficient  B  were  set  equal  to  zero. 

In  presenting  numericed.  results,  solid  curves  have  been  used  for 
plotting  results  pertaining  to  the  strain  hairdening  materiaJ.,  whereas 
the  corresponding  plots  for  the  perfect  plastic  medium  are  indicated 
by  broken  curves. 
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7-2  Wumerlcal  Results  ajid  Discussion  for  Plate  Problems  auad  for  the 
HBLlf -Space  Problem 

The  damage  inciirred  in  a  half -space  or  plate  by  a  heat  pulse  may 
be  characterized  by  the  magnitudes  of  residual  stresses  and  the  dimen¬ 
sion  of  the  total  plastically  deformed  region.  As  far  as  the  residvial 
stresses  are  concerned,  the  location  and  the  largest  magnitude  of  the 
residual  stress  are  of  great  importance,  and  will  be  considered  first. 
We  recall  frem  the  dlsciission  presented  in  the  preceding  section  that 
the  maximum  residual  stress  occurs  at  the  free  boundaries  of  both  a 
half -space  and  a  plate.  Moreover,  the  largest  residual  stresses  in  a 
half -space  and  a  plate  are  eq\ial  if  the  associated  maximum  boundary 

temperatures  and  T  -  are  the  same.  The  maximum  residued.  stress - 
m  mL 

es  were  calculated  for  different  values  of  T  or  T  ^  for  both  a 

m  mL 

strain  hardening  material  having  a  tenqperature  dependent  yield  stress 
and  a  perfectly  plastic  material  having  a  constant  yield  stress.  This 
calculation  was  based  on  the  formulas  (If. 5.12),  (7.1.5);  (5.5«8),  or 
(4.6.11),  (7.1.4),  (5.5.18),  depending  on  whether  plastic  flow  occurs 
during  unloading  or  not.  The  results  of  the  calculation  for  values  of 
^m  "^mL  from  175°P  to  700°F  are  plotted  in  Fig.  7-1.  The 

critical  temperatures  T^^  and  T^  for  which  plastic  flow  Impends 
in  a  half -space  and  a  plate,  respectively,  were  calculated  from  (4.4.6), 
(6.5.7)  and  the  material  properties  given  in  (7. 1.5).  The  temperature 
specifying  incipient  plastic  flow  in  unloading  was  caloxlated  either 
from  (4.6.3)  or  from  (6.5.7).  It  was  found  that 


’^1  -  ‘ 


575.06°p 


(7.2.1) 
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for  the  strain  hardening  materleuL,  and 

■'i  ■  V  '  ’«  ■  ’k  ■ 


(7.2.2) 


for  the  perfectly  plastic  material. 

We  note  from  Fig.  7-1  that  the  mEucimum  residual  stress  is  linear¬ 
ly  related  to  the  maximum  temperature  or  T^.  Actually,  the  lin¬ 
earity  between  the  mEiximum  residual  stress  and  the  maximum  tesi^erature 
can  be  readily  observed  on  inspection  of  (4.5*12),  (4.6.11)  or  (7.1. 5)> 
(7.1.4).  In  view  of  (7-2.1)  and  (7.2.2),  we  further  note  from  Fig.  7-1 
that  the  maxlmm  residual  stress  for  the  strain  hardening  material,  in¬ 
creases  more  rapidly  with  the  maximvim  temperature  for  values  of 
less  than  than  for  values  of  greater  than  Here  !I1^ 

corresponds  to  the  maximum  temperature  for  which  incipient  plastic  flow 
dxirlng  unloading  Impends.  If  the  maximum  tenperatxire  is  greater 
than  the  critical  vadue  the  material  is  in  the  plastic  state  in¬ 

stead  of  in  the  elastic  state  as  being  characterized  by  vailues  of 
less  then  the  critical  value  This  accounts  for  the  fact  that  the 

curves  relating  the  maximum  residual  stress  end  maximum  tenQ>erature 


change  their  sieves  at  The  foregoing  remarks  make  it  clear 

that  the  maximum  residuel  stress  for  the  perfectly  plastic  material 


would  remain  constant  for  T  greater  than  the  critical  value 

m 


(Fig.  7-1). 


The  dimension  the  total  plastically  deformed  region  is  also  a 

characteristic  measure  of  the  (iamage  incurred  in  a  helf -space  or  plate 
by  a  heat  pulse.  Therefore,  we  n(3W  proceed  to  calculate  the  position 
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of  the  boundary  of  the  total  plastically  deformed  region.  The 
positions  R^  and  R^j^  of  the  total  plastically  defomed  and  steady 
state  plastic  regions  in  a  half -space,  suid  corresponding  to  different 
values  of  T^^  may  be  calculated  from  (4.6.4)  suid  (7.1.5)*  In 

order  to  deteimine  R^  and  R^  in  a  plate,  the  steps  suggested  in 
section  6-3  were  followed.  The  dimensionless  values  of  R^^  for  y 
equal  to  0.5  are  plotted  in  Fig.  7-2  both  for  a  half -space  and  a  plate, 

7  being  a  parameter  introduced  in  the  ten^perature  functions  (6.2.1) 
and  (6.2.2)  to  determine  the  shape  of  the  heat  pxxlse.  For  a  plate, 

7L  represents  the  distance  between  the  source  plane  and  the  free  bound¬ 
ary  of  the  plate  (Fig.  6-l).  Similar  plots  showing  the  variation  of 

both  the  position  R,  eind  with  mEocimum  boundary  temperature  T  » 

X  XX  2QL 

for  the  individxial  plates  for  values  of  y  equal  to  0*50  and  1.00  are 
presented  in  Figs.  7-3  and  7-4. 

In  the  case  of  the  plate  having  a  fixed  and  insulated  boundary, 
two  critical  temperatures  T^^  suid  T^^  ^ere  found  to  exist.  Speci¬ 
fically,  T^^  is  the  maximum  boundary  tenqperature  for  which  the  botmd- 
ary  of  the  total  plastically  defoimed  region  coincides  exactly  with  the 
insTilated  boundary  of  the  plate,  whereas  T^g  is  the  maximian  boundary 
tenqieratxzre  for  which  the  boundary  of  the  steady  state  plastic  region 
coincides  with  the  insvilated  boundary  of  the  plate.  The  values  of 

T  ,  and  T  „  are  deteimined  by  setting  each  of  the  quantities  T  „  / 
cx.  Cc 

^cl  '^02  '^mc/’mL*  '^mR^  ^ML 

ten^eratures  for  which  plastic  flow  ijaqiendB  during  loading  and  imload- 

ing,  respectively,  and  T^^^^  is  the  maximum  tesiq)erature  attainable  at 

the  Insulated  boundary. 
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200  300  400  500  600  700 

Temperoture  or  T^l_ 

FIG  7-2  THE  DIMENSION  OF.  THE  TOTAL  PLASTICALLY  DEFORMED  REGION  WITH 
TEMPERATURE  Tm  or  TmL 


FIG  7-3  VARIATION  OF  THE  DIMENSIONS  OF  STEADY  STATE  PLASTIC  REGION  AND  TOTAL 

PLASTICALLY  DEFORMED  REGION  WITH  MAXIMUM  TEMPERATURE  OF  FREE  BOUNDARY 
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FIG.  7-4  VARIATION  OF  THE  DIMENSIONS  OF  THE  STEADY  STATE  PLASTIC  REGION  AN 
THE  TOTAL  PLASTICALLY  DEFORMED  REGION  WITH  MAXIMUM  FREE  BOUNDARY 
TEMPERATURE. 
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For  the  numerical  values 

r  .  0.50  ,  .  0.6671.  , 

ve  obtain 

=  266.7066®?  ,  =  56i.9718°F  ,  (7-2.5) 

for  the  strain  hardening  medium  and 

T^l  =  3^.6209®F  ,  T^g  =  689.21*.i8®P  ,  (7-2.4) 

for  the  perfectly  plastic  medium. 

The  next  calculation  concerns  the  stress  and  deformation  of  the 
half -space;  the  transient  positions  of  the  interfaces  pg  and  p^ 

were  calculated  from  (4.4. 5)>  (4.5-l)>  (4.6.26)  and  (5-2.27).  The  re¬ 
sults  are  given  in  Fig.  7-5»  the  maximum  temperature  being  used  to 
label  the  c\irves. 

Next,  the  equations  (4.4.2),  (4.5-2),  (4.6.2a)  were  used  to  deter¬ 
mine  the  instants  t, ,  t  _  and  t  for  the  values  300°F,  400°P, 

J.  n®^  pa 

500®F,  600°P  end  700°F  of  the  meucimum  boundary  temperature  T^^^.  The 
dimensionless  results  of  this  conqputatlon  together  with  end 
for  the  above  mentioned  temperatures  are  presented  in  Table  7”1- 

It  was  shown  in  chapter  4  that,  if  the  maximum  boundary  tenq>erature 
Tjjj  is  greater  than  the  value  !I!^  given  by  (4.6.5),  and  if  the  instant 
t  -  is  greater  than  the  instant  t  ,  then  four  regions  exist  simultan- 

Buu  PCL 

eously  between  the  instants  t  and  t  _  in  the  half -space.  This  may 

j)ci  nui^ 
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FIG  7-5  -POSITIONS  OF  DIFFERENT  INTERFACES  AT  DIFFERENT  INSTANTS 
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R^/a  and.  R^.j/a  denote  the  dimensionless  positions  of  the  boundary  of  the  totaJ.  plaatloally  defoned 
and  steady  state  plastic  regions,  respectively,  n^  and  n^  are  defined  as  the  dimensionless  tine  Instants 
at  vhleh  plastic  flew  during  loading  impends  at  z  -  a  and^  x  ■  R.  ;  n  is  defined  ais  the  dimensionless 

X  pft 

tlm  instant  at  vhich  incipient  plastic  floir  daring  imloadlng  occnars  at  x  «  a. 
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be  verified  in  Pig.  7-5  for  values  of  T  equal  to  600  F  juid  700  F  for 

m 

the  strain  hardening  material,  suid  for  equal  to  700  P  for  the 
perfectly  plastic  material. 

We  observe  from  Fig.  7-5>  or  from  (i4-.4.7)>  that  the  same  maximum 
boundeury  temperature  will  affect  plastically  a  Isurger  region  if  the 
yield  stress  is  temperature  dependent.  The  same  observation  is  valid, 
of  course,  also  for  plates. 

Making  vise  of  Fig.  7-5#  and  letting  =  4CX)  F,  stresses  and 
deformations  were  calculated  from  (4.5.9),  (4.3»8),  (4.5«15)  for  the 
elastic  region,  from  (4.4.26),  (4.4.27),  (4.4.28),  (5.5.1),  (5.3.2), 

(5.3.3)  for  the  plastic  regions,  from  (4.5.5),  (4.5*6),  (4.5*7), 

(5.3.4) ,  (5.3.5),  (5.3.6)  and  (5.3.7)  for  the  elastic  region  of  tmload- 
ing  in  both  the  strain  hardening  amd  perfectly  plastic  media.  The 
dimensionless  results  of  the  calculation  are  presented  in  Figs.  7-6  to 
7-11.  As  time  goes  on,  the  transient  stresses  and  deformations  tend  to 
their  steady  state  values  which  were  computed  from  (4.5*12),  (4.5.15), 
(4.5.14),  (5.3.8),  (5.5.9),  (5.5.10)  for  the  steady  state  elastic  region, 
and  from  (4.6.11),  (4.6.12),  (4.6.15),  (5.3.8),  (5-3.9),  (5.3.10)  for 
the  steady  state  plastic  region.  The  dimensionless  resvlts  of  the 
coogjutation  for  different  values  of  the  naximum  bovmdary  tenqjerature 

T  for  the  residual  stress,  strain  euid  displacement  are  shown  in 

SI 

Figs.  7-12,  7-15,  and  7-14,  respectively. 

Fig.  7-6  shows  that  the  largest  ccnq)ressivB  stress  Induced  in  the 
plastic  region  is  greater  l^t  the  perfectly  plastic  medium  than  in  the 
strsdn  hardening  medium  having  a  temperature  dependent  yield  stress. 

The  same  is  not  true,  however,  for  residual  tensile  stress.  This  is 
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FIG.  7-6  VARIATION  OF  STRESS  WITH  TIME  AT  DIFFERENT  POSITIONS  IN  THE  HALF  SPACE 
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FIG  7-7  VARIATION  OF  STRAIN  WITH  TIME  AT  DIFFERENT  POSITIONS  IN  THE 
HALF  SPACE. 
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FIG  7-e  VARIATION  OF  DISPLACEMENT  WITH  TIME  AT  DIFFERENT  POSITIONS  IN  THE  HALF  SPACE 
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FIG  7-9  STRESS  DISTRIBUTION  AT  DIFFERENT  TIME  INSTANTS  IN  THE  HALF  SPACE 
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FIG  7-10  STRAIN  DISTRIBUTION  AT  DIFFERENT  INSTANTS  IN  THE  HALF  SPACE 
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FIG  7-13  RESIDUAL  STRAIN  DISTRIBUTION  IN  THE  HALF  SPACE 
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due  to  the  fact  that  the  txanslent  stresses  In  hoth  the  strain  hardening 
and  perfectly  plastic  media,  after  falling  back  from  their  msoctmum  value 
In  ccmpresslon  to  zero,  begin  to  increase  In  the  opposite  direction;  the 
transient  tensile  stress  In  the  strain  hardening  medium  is  at  any 
Instant  hl£^er  than  that  in  the  perfectly  plastic  medim  (Fig.  7 "6)  and 
eventually  assmes  a  steady  state  vEQ.ue  also  higher  than  the  one  assumed 
by  the  perfectly  plastic  medium.  The  transient  emd  residual  strains  and 
displacements  are,  hovever,  greater  In  the  strain  hardening  medium,  as 
may  be  seen  from  Figs.  1-T,  7-8,  7-10,  7-ll>  7-13,  and  7-1^.  Another 
point  worth  noting  is  that  in  both  the  transient  and  steady  state  solu¬ 
tions  (Figs.  7-9  to  7-1^)  the  stresses,  strains  emd  displacements  are 
always  greatest  at  the  boundary.  We  further  observe  from  Figs.  7-12  and 
7-13  that  the  residual  elastic  stress  smd  strain  increase  much  more 
rapidly  with  the  maximum  boundary  tenqperatiire  than  their  plastic 
counterparts  in  the  total  plastically  deformed  region;  the  position 
of  the  steady  state  elastic  plastic  interfaces  can  be  read  off  from 
these  figures  by  locating  the  planes  at  which  there  is  an  abrupt  change 
in  the  slope  of  the  curves.  For  =  400  F,  the  maximum  transient  and 
residual  stresses,  strains  emd  displacement  were  found  to  be  as  follows: 

T y/y^  m  -  1.5488  ,  at  n  =  n^  =  9.2485  , 

€ Tjjj  =  2.5840  ,  at  n  =  1.00  ,  (7-2.5) 


Sj/o  -  0.8a95  .  Oj/ao  \  -  4.0147 
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for  the  strain  hardening  medium  and 


Ty/y^  -  -  2.CXX) 

,  at  n  =  n^^  ®  0.3050  , 

\  “  2.2860 

,  at  n  =■  1.0000 

(7.2.6) 

Tjj/y^  =  1.4800 

,  T,  .  o.sw  ,  Ug/m  . 

-  5.7950 

for  the  perfectly  plastic  medium. 

Recalling  the  discussion  presented  in  section  7-l>  we  conclude  that 
the  values  given  by  (7.2.3)>  (7*2.4)  for  the  nmimum  transient  and 
residual  stresses  and  strains  also  represent  the  corresponding  quanti¬ 
ties  induced  in  a  plate  by  a  heat  pulse  of  the  same  maximiBn  amplitude, 
i.e.,  for  =  400  F. 

The  transient  displacements  are  shown  in  Fig.  7-8;  the  fact  that 
even  in  the  elastic  region  the  residml  displacement  does  not  vanish 
is  due  to  the  regularity  condition'^’used.  Therefore,  in  presenting  the 
permanent  deformations  in  the  half -space,  the  quantity 

[Ug/aa  +  (l+l/^/l-2/ ) 

was  plotted  against  position  x/a  (Fig.  7-14).  The  resulting  curves 
are  found  to  decrease  from  a  naximum  value  at  the  boundary  of  the  hed.f- 
space  to  zero  at  the  boundary  of  the  total  plastically  deformed  region. 
7-5  Further  Discussion 

In  the  preceding  sections  of  this  chapter  we  detemined  the 
dimensions  of  the  total  plastically  deformed  emd  steedy  state  plastic 
regions  and  slIso  obtained  expressions  for  the  resldiiEd  stresses  and 


I  cf.  see  equation  (4.3*6) 
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strains  caused  by  a  uniformly  applied  heat  pulse  to  the  boundary  of  a 
half -space  and  a  plate.  We  are  nov  confronted  with  the  problem  of 
predicting  the  response  of  the  medium  to  an  additional  heat  pulse 
applied  after  conditions  produced  by  the  first  pulse  have  reached  a 
steady  state.  In  particular^  we  wculd  like  to  find  the  amplitudes  of 
the  second  heat  pulse  for  which  plastic  flow  will  recur. 

Let  us  consider  the  application  of  a  second  heat  pulse  in  seme 
detail,  assuming  that  steady  state  condition  has  been  reached  for  the 
first  pulse.  Hie  relation  between  maximum  shearing  stress  q  emd 
temperatvtre  T  Is  given  by 

.  0$- 
1  =  ^ 

where  the  maximum  shearing  stress  q  will  Increase  from  Its  residual 
value  ^  as  the  temperature  T  Increases  from  Zero.  Thus  integration 
of  (7.^.1)  yields  the  relation 

1  =  iOTzp  T  +  qj,  .  (7.4.2) 

Recalling, 

q  =  -  =  "  -#  »  (4.3.12)  R 

we  may  write  (7.4.2)  as 

aE  T 

“  2(l-Z/j^)  ^  ■  2 


<1 


(7.4.3) 


'rtiere  is  the  resldtial  stress  for  which  the  distribution  in  the 
half -space  is  shown  in  Fig.  7-12.  Plastic  flow  due  to  biaxial  ccw- 
presslon  will  resume  wherever  the  maximum  shearing  stress  q,  given  by 
(7*^«3)  beccaies  equal  to  the  current  yield  stress  y^  in  shear,  where 

=  Fo  -  j—sup  p"  (7.4.4) 

The  relation  (7.4.4)  is  actually  the  relation  (f)  in  section  4-42. 


Fig.  7-15  Shearing  Stress,  Strain  and  Current  Yield  Stress 
From  Fig.  7-15  we  observe  that 
P"  '  Pj,  +  p’  , 

Pp  teing  the  principal  residual  shearing  strain. 


If  we  recall  that 


(4.3.12)  R 


then  the  expression  for  the  current  yield  stress  becomes 

=  Yo  -  BT  +  (2up  -f  +  8  -|  )  .  (7-4.5) 

Replacing  q  in  (7.^.3)  by  given  by  (7.4. 5)>  we  find  the  terapera- 
tiire  above  which  plastic  flow  will  again  occur  as  the  result  of 

the  second  heat  pulse: 


■■lx 


A 


B  +  aE^/2(l-  y^  1-0^  1-8 


i]  .  (7.4.6) 


As  may  be  readily  seen  from  (7.4.6),  the  ten5>erature  is  not 

constant  but  varies  with  the  residual  stress  and  strain  in  the  total 
plastically  deformed  region.  In  the  region  that  has  never  been  plasti¬ 
cally  deformed,  the  residual  stress  and  strain  are  zero,  and, 
in  this  case,  (7-4.6)  reduces  to  the  form 


^Ix'  “  B  +  0!E^/2(l-^^)  (7.4.7) 

Identical  with  (4.4.6).  Therefore,  as  long  as  the  medium  remains  in  the 
virgin  state,  the  temperat\Jre  above  which  plastic  flow  is  incipient 
remains  equal  to  T^,  regardless  of  the  number  of  pulses  to  which  the 
medium  may  have  been  subjected. 

With  the  valves  of  the  residual  stress  and  strain  given  in 
Figs.  7 "12  and  7 ”13  for  the  half -space,  and  the  values  of  the  maximum 


residual  stress  and  strain  obtainable  from  Figs.  T-l  and  7-15  for  tne 
plate,  the  temperature  for  the  strain  hardening  meditan,  and 

corresponding  to  different  values  of  T  or  T  ,,  was  calculated,  the 

la  ULL 

results  being  summarized  in  Table  7-2. 


Table  7-2  Ten^ieratures  Corresponding  to  Recurrence  of  Plastic 


Flow  at  the  Heated  Boundary  of  Half -Space  or  Plate 


300°F 

4CX)°F 

500°F 

6qo°f 

700°p 

178 

178 

178 

178 

178 

Ti,'  (°F) 

509.9692 

425.0640 

The  values  of  Tj^^  for  other  positions  in  the  medi\im  can  also  be 
determined  readily  by  direct  substitution  of  the  corresponding  values  of 
T  and  €  into  (7.4.6).  Ccraparing  (4.4.6)  with  (7.4.6),  we  observe 

It  K 

tliat  as  a  resxilt  of  the  first  heat  pulse,  the  temperature  corresponding 
to  the  incipient  yielding  in  the  plastically  deformed  region  has  been 
raised  by  a  factor  of 


[1  +  1/2  (.1  4- 

'  y  l-s  l-s  y 


Vz  further  note  from  the  results  presented  in  Table  7-2  that  if  the 

maximum  amplitude  T  or  T  ,  of  the  first  heat  pulse  remains  within 
m  ml) 

certain  limit,  no  further  plastic  flew  will  be  Induced  by  the  second 

heat  pulse  having  the  same  amplitude  as  the  first  one.  '-Che  case  where 

T  =  500  F  may  serve  to  illustrate  this  point  (Table  7-2).  On  the  other 
m 

hand,  if  the  nBxtmum  amplitude  of  the  first  pulse  exceeds  this  limit, 
and  the  maximum  amplitude  of  the  second  pulse  is  equal  to  that  of  the 
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first  ptilscj  further  plastic  flov  will  be  Induced,  and  plastic  strains 
will  tend  to  accumulate.  This  contention  is  borne  out  by  an  inspection 
of  the  cases  where  or  T  -._  is  eqiusil  to  400°F,  500°P,  600°P  cuid 

in  nUi 

700°p  (Table  7-2). 

Por  a  more  detailed  determination  of  the  plastic  flew  during  a 
second  heat  pulse  a  coc^lete  anGilysls,  ccmparable  to  that  presented  for 
the  first  heat  pulse,  is  needed.  Although  we  have  not  pursued  this 
analysis,  the  procedures  outlined  in  this  work  are  applicable  and 
sufficient  for  the  study  of  any  subsequent  heat  ptilses. 
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implied  Mechanics. 

4.11  C.  Huang  "Thermal  Stress  in  an  Elastic  Work  Hardening  Sphere" 
Journal  of  APPH^J  Mechanics,  Vol.  27,  Series  E,  No.  4, 

pp  629-634,  i960. 

4.12  H.  S.  Carslaw  and  J.  C.  Jaeger  "Conduction  of  Heat  in  Solids" 

2nd  Edition,  Oxford  Press,  p  259>  1959* 

6.1  Same  as  4.4. 

6.2  Same  as  4.6. 

6.3  Same  as  4.7* 

6.4  Same  as  4.3* 

6.3  Same  as  4.8. 
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State  Plastic  Re^on  in  Plates.  (Figs.  7-5t  7-4). 


Fig.  A-1  For  the  Plate  ¥ith  the  Boundary  at  x  >  L  Subjected  to  a  Heat 
Pulse  and  With  the  Other  Boundary  at  x  >  o  Maintained  at 
Zero  Teo^erature. 


Fig.  A-2  For  the  Plate  With  the  Boundary  at  x  »  L  Subjected  to  a  Heat 
Pulse  and  With  the  Other  Boundary  at  x  >  o  Insxilated. 
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FIG  A-l  VARIATION  OF  MAXIMUM  TEMPERATURE  AND  CORRESPONDING  INSTANTS  ALONG  PLATE 


Numbers  in  parenthesis  denote  values  of  y  in  16-2-2)  and  (6- 2-5) 
Tmxond  Dmx  curves  show  maximum  dimensionless  temperature  ond 
corresponding  time  instants  along  plate 
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FIG.  A-2  VARIATION  OF  MAXIMUM  TEMPERATURE  AND  CORRESPONDING  INSTANTS  ALONG  PLATE 
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APPENDIX  B 


Derivation  of  Plastic  Stress  Strain  Relations  s  emd  s  Given  By 

— —  I  ■  -  I  I  XX'-  ' ■  yy  ~  ■  -  ■ 

(a)  and  (t>)  in  (4.4.14) 

The  relation  (4.4.13)  in  principal  stress  space  nay  be  written  u 


s^  ■  E  * 

r-1 


(7) 


(4.4.15)  E 


(dF^’'^/&Vjjj  Vjj  +  aF^’'^/aT  T) 


the  yield  functions  are  given  by 


/, \  V  -  V 

i.(l)  „  JS _ L  , 


2 - (»'2-®+  - 


(4.4.8)  R 


.(2)  _  JLlZ^ 


2— . 


(4.4.9)  R 


/,  s  V  -  V  T 

,(1)  .  _x _ JL  .  _  _22: 

2  2* 


J2)  \  ■  ""z  ""zz 

f '  »  2 - -  “  "2"  ^ 


(4.4.10)  R 


Renal  li  ng  that 


T 

XX 


0 


T  =  T 

yy  zz 


and 


ij  ij 


- 1/^  •> 


ij  kk 


we  have 


By  (k.k.8),  (4.1»-.9)y  (4.4.10)  and  (a),  we  obtain 
^^^^/^XX  “  V2  ,  =  -  1/2  , 

(b) 

=  1/2  ,  -  1/2  ,  =  0  , 

X  y  z 

euid 

=  1/2  ,  =  - 1/2  , 

(c) 

=  1/2  ,  ap^^^/^^  =  0,  » - 1/2  . 

X  y  z 

By  (4.4.8),  (4.4.9)  anfi.  (4.4.10)  it  was  found  that 

SH/ai  -  (lA-s)/!/?  »^p  > 
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By  («•)>  (B)  and  (d),  the  following  resxilts  are  obtained: 

1/2 


,  «  1/2  1/2 
[1/2  1/2  +  (-  l/2)(-  1/2)  +0  0]  =  (1/2)  ,  (e) 


Vj^  +  T) 

=  +  bT^^^/dv^  Vy  +  v^  +  aF^^^/&T  T) 

=  [1/2  (-  2/3  Tyy)  +  (-  l/2)(l/3  Tyy)  +  0  (1/5  Tyy)  +  HP] 

=  -  1/2  +  HE  .  (f) 

Similarly,  It  was  foxind  that 

(st<=>/aT^  -  (1/2)^'^^  (g) 

wd 

(3f^^^/^j^  Vj^  +  bF^^^/bl  T) 

"  -  1/2  Tyy  +  HP  (h) 


f 
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Substituting  (b),  (c),  (d),  (e),  (f)^  (g),  and  (h)  in  (1^.4.13)> 
we  finally  obtain 


8  =8 
XX  XX 


(X)  ,  ,  (2) 


XX 


•  • 


(  -  1/2  Tyy  +  OT) 


(a)  in  (4.4.14) 


.  W  .  ,  -s 

yy  yy  zz  zz 


(b)  in  (4.4.14) 
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APPENDIX  C 


Redaction  of  the  Strain  Hardening  Solution  to  Perfectly  Plastic  Solution 

It  will  be  shown  In  what  follows  that  (3 •2.6)  Is  reducible  from 
(4.4.16). 


Ve  recall  that 


and 


4(1-22/^) 

~ 


B  +  30]  T  , 


€  =  N  T  , 

X  ' 


where 


jj  ^  ^  B„il,:g.laO,,.]  .  M  .  -ii] 


2(x. 


E 


M 


2B  E^  (1-s)  -  2  Ej^ 

Ej^  (1-s)  +  2(1- y^)  u  /ITT 


(5.2.6)R 

(4.4.16)R 


(4.4.17)R 


We  now  define 

4(1-2^) 

"pp-t - (a) 


« 


If  In  (4.4.17)  la  set  eqial  to  zero,  N  should  be  reducible  to 

N  .  N  may  be  rewritten  as 
PP 


-=-iM  +/T  (1-s) 
®1 


2B  -  M 
2  h. 


(b) 


5-0- 
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Substituting  M  given  by  (1(-.4.17)  in  the  quantity  (28-*!)/2Up, 
we  eventually  have 


2B  a 

E^(1-8)  +  2(1-2^^)  hp/175 


(c) 


If  Up  is  set  equal  to  zero,  corresponding  to  a  perfectly  plastic 
material,  the  foUcving  results  are  obtained. 


M  a  2B  ,  s  =  hpAg  =  0 
2B.M  ^  2B(l-i/^)^E,a 


2  h. 


=  /V5 


=1 


(d) 


Now,  if  (d)  is  substituted  back  in  (b),  we  obtain 


2B(1-2/) 

K  =  a  -  2B  +  /T  /I73  [ - g—i-  +  O] 


hp  Kl-U) 

=  o  -  B  - - ^  B  +  2  a 


k(i-  2t/^) 
E 


B  +  5  a]  =  N, 


pp 
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